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I. 



INTRODUCTION 



The advent of the digital computer as standard 
eguipment cn hoard virtually all modern Naval ships has 
opened the field of Digital Computer Control in almost all 
aspects cf ship life. The computer has been a viable asset 
in fire ccEtrol systems for years and has been used 
extensively for aids to ship maneuvering control in the form 
of HIDS (Naval Tactical Data System) readouts. The 
declining costs cf general and special purpose computers has 
made their inception as a manpower replacement or 
augmentation a reality. Their high speed and accuracy can 
make them perform functions with much greater safety than 
previously attainable with time proven (and sometimes time 
weary) "seaman's eye." 

This then is the basis for this thesis; a study of the 
types of maneuvers that can be handled more accurately and 
safely than presently being accomplished. 

One area of study is the total Replenishment At Sea 
(RAS) problem including both course and speed control for 
the approach and alongside phases. This situation has 
always been one of extreme danger due to the collision 
potentials involved. However, other dangers are involved in 
the on deck evolutions when the ships are not kept at a 
fairly constant distance. Sudden violent maneuvers may 
cause the stress on the lines between ships to increase 
enough to cause the lines to part. The reality of this 
danger is readily apparent to anyone who has ever seen a 
marila line or steel cable part or a kingpost shackle break 
or a kingpost suddenly bend under these extreme stresses. A 
system which will minimize these potential dangers is well 
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wcrth investigation. 



Of course with a digital computer, the algorithm for 
BAS can he modified or replaced to enable its use as a 
maneuvering control device for other situations such as 
formation steaming or single ship navigation transit 
control. 
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II. MODELING 



A. MARINES EINAMICS 

In the conception of this thesis, realistic models of 
modern destroyer hull configurations were sought. This 
search proved fruitless. The hydrodynamic coefficients for 
present cay destroyers are not currently available. 
However, seme naval and civilian research is presently being 
conducted to obtain these coefficients. 

A complete set of these coefficients is necessary for 
any maneuvering control system design. A hull configuration 
which has been under continual study with well defined and 
verified hydrodynamic coefficients was chosenC 1 ]., This hull 
is commonly referred to as the ,, Mariner ,, hull. 

The development of the equations of motion in six .and 
three degrees of freedom have been well documented [ 2 ] . The 
model used for this thesis is the equations of motion in 
three degrees of freedom linearized with second order and 
higher terms eliminated. These equations are characterized 
by dependency on small perturbations about a specific 
operating point. The maneuvers experienced in the following 
chapters dc not entirely meet this criterion. The 
inadequacy and shortcomings of this model are of little 
ccnseguence because relevant hydrodynamic coefficients are 
not available, and the methods presented can be applied to 
any ship type. 

The development of the model is readily available to 
the interested reader [ 3 ] ; only a summary of the equations 
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and tfaeix corresponding hydrodynamic coefficients 
presented here. 

The eguations of motion used are as follows: 

(X- -it) U+X (0-u ) +X d=0 

u u 0 d 

(Y.-n) Y+I V+(Y -m)R+Y.R+Y d=0 
v v r r d 

(N.-I ) E+N R+N.V+N V+N d=0 

r z r v v d 



The direction and sense of the terms in the 
equations are shown in figure II- 1. Letting: 



a =m-Y. 
11 v 

b =-Y 

11 v 

c =C 

11 

a =-Y. 

21 r 

b =m- Y 

21 r 

c =C 
21 

a =-N. 

12 v 

t =- N 
12 v 

c =C 
12 

a =1 -S. 

22 z r 

b =-N 
22 r 

c =0 
22 

a =m-X. 
33 u 



are 



above 
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X 



Figure II-l 

Direction and Sense of Hydrodynamic Terms 
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b =-X 
33 u 

c =0 
33 



Setting 

V=A 

V=B 

U=C 

The equations cf motion can be written a 



n •• 

a A+b A+c A + a B+b B+c B=IF1 

11 11 11 21 21 21 

« it 

a A+b i+c A+a B+b B+c B=IF2 

12 12 12 22 22 22 



a C+b C+c C=IF3 
33 33 33 



cr : 



n n 

a A+a E=I 1 
11 21 

ii ti 

a A+a E=I2 
12 22 



a C=I3 

7 7 



where: 



IF 1=-Y «d 
d 



IF2= N «d 
d 
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IF 3 = 




X 

u 



dt 



and : 



I1=-t 1-c A-b B-c B + IF1 

11 11 21 21 



I2=-1 l-c A-b B-c B+IF2 
12 12 22 22 



l3=-t C-c C+IF3 
33 33 



By solving this system of equations, the following 
relationships are established: 



A= (a a I1-a a I2)/DEI 
22 33 33 21 



E= (a a I2-a a II) /DEI 
11 33 33 12 



it 

C=I3/a 

3 



where; 



DEL=a (a a -a a ) 
33 11 22 12 21 



which yields the solution: 



v=iUv 

0 



U=B = « 

0 




II 

A dt 



B 



d t= 





U=C=u 

0 




ii 

C dt 



ii 

Edt ]dt 



The space coordinate system is defined as fellows: 



23 



I=0«SIN(| + V*C0SP 



x=o«cosp-T?*smp 



where : 




i dt 



i dt 



Table II- 1 ? summarizes the symbols and nomenclature used 
in the foregoing abbreviated solution of motion in three 
degrees of freedom. The applicable characteristics of the 
Mariner hull are presented in table II-2 with the 
corresponding nondimensionalized hydrodynamic coefficients 
and the DSL computer program variable names delineated in 
table II-3. 



Computer Program #1 is the basic DSL program that was 
developed from these eguations of motion. This program uses 
two ships to illustrate the turning characteristics of the 
Mariner bull for various rudder commands. Figure II-2 shows 
the difference between a step model and a ramp model rudder 
in a geographic plot.. Figure II-3 shows the corresponding 
difference in yaw. 



24 



TABLE II-1 



Sym bol 



X. 

u 



0 

a 



Y. 

v 



Y 

v 



Y. 

r 



Y 

r 



Y 

d 



N . 

v 



N 

v 



SYMBOLS AND NOMENCLATURE 



Def in ition 



derivative of longitudinal force with 
respect to longitudinal acceleration 0 

derivative of longitudinal force with 
respect to longitudinal velocity U 

derivative of lateral force with respect to 
transverse acceleration V 



derivative of lateral force with respect to 
transverse velocity V 



derivative of lateral force with respect 
tc angular acceleration R 



derivative of lateral force with respect 
to angular velocity R 



derivative of lateral force with respect 
to rudder angle d 



derivative of yaw moment with respect to 
transverse acceleration V 



derivative of yaw moment with respect 
tc transverse velocity V 
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TABLE II-1 (cont.) 





SYMBOLS AND NOMENCLATURE 


Symbol 


Definition 


N. 

r 


derivative of yaw moment with respect to 
angular acceleration R 


N 

C 


derivative of yaw moment with respect tc 
angular velocity R 


N 

d 


derivative of yaw moment with respect to 
rudder angle d 


• 

B 


yaw angle acceleration 


B 


yaw angle velocity 


u 

0 


initial velocity of origin cf body axes 
relative to fluid 


• 

7 


transverse acceleration of ship axes 
relative to fluid 


7 


transverse velocity of orian of ship axes 
relative to fluid 


X 


hydrodynamic longitudinal fcrce 


Y 


hydrodynamic lateral force 


• 

0 


longitudinal acceleration cf ship axes 
relative to fluid 
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TABLE II- 1 (cont.) 



SYMBOLS AND NOMENCLATURE 



Sym bol Definition 

U longitudinal velocity of ship axes relative 

to fluid 

^ yaw angle 

At actual time 

T nondimensionalized time 

x longitudinal distance that the ship center 

g 

of gravity is forward of the ships axes 

u longitudinal velocity of ship axes relative 

tc fluid (operating point) 



TABLE II-2 

CB ABAC T ERISTICS OF MARINER-TYPE 


STUDY SHIP 


Length, ft. 


527.8 


Beam, ft 


76.0 


Draft, ft 


29.75 


Displacement, tons 


16,800. 


Block coefficient, C 


0.6 
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TABLE II-3 



NCNDIME NSICN AL HYDRODYNAMIC COEFFICIENTS 





Computer Proqram 


Non dimensional 


:oef icient 


Variable Name 


Value 


(X. -IB) 
U 


MXOD 


-0.0085 


X 

u 


XU 


-0.0012 


Y 

V 


YV 


-0.01243 


(Y . -in) 

V 


MYVD 


-0.015 


(Y -mu ) 
r 1 


MYR 


-0.0051 


(Y.-mx ) 

r g 


YRD 


-0.00027 


Y 

a 


YDELR 


+0.0027 


N 

V 


NV 


-0.00351 


N. 

V 


NVD 


-0.000197 


(N -mx u ) 
r g 1 


NR 


-0.00227 


(N.-I ) 
r z 


IZNRD 


-0.00068 


N 

d 


NDELR 


-0.00126 


X 

□ 


— 


-0.0000462 


Y 

□ 





-0.0000052 


N 

n 


— - 


+0.0000026 


X 


XDELR 


0.0 



d 



NOTE: X =C.O 

g 
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TABLE II-3 (cont.) 



NCNEIME KSIONAL HYDHOD YNAMIC COEFFICIENTS 

Values tased on the following operating point: 
a =1.0 (15 Kts ) 

i=o.o 

v=o.o 
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Rudder Step and Ramp Model Geographic 



2 
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Rudder Step and Ramp Model Yaw 
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E. 



BUDDIE BESP0N5E 



The previous section indicates a marked difference in 
behavior between step and ramp rudder models. This prompted 
an investigation into realistic rudder modeling which would 
fulfill the reguirements of limit stops and maximum rudder 
rate. 



NSBEC[ 4 ] has modeled the rudder of the DD-931 Class 
Destroyer. The basics of this model are presented in the 
block diagram of figure II-4. 

The first limiter models the rudder steps which for the 
Mariner are ±30 degrees. The second limiter models the 
proportional band of a variable-displacement pump by 
limiting its maximum percent stroke. The limits for this 
nonlinear element have been found to be ±7 degrees. 

The transfer function (K /s) accepts an input error 

9 

signal of up to 7 degrees, converts it to a rudder rate, and 
integrates the rudder rate to obtain rudder angle. letting: 

d » Maximum rudder rate (2.0 degrees/sec) 
m 

d' = Maximum error input (7.0 degrees) 
emax 

The system gain can be defined as: 

K = d /d 
g m emax 

= 2.0/7. 0 

= 0.265714 /sec 
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Tc convert this model to the required nondimensionalized 
form, the following manipulation is required: 



K ' = K •L/u 
g g 1 

= 0.285714*527.8/(15*1.689) 
= 5.95224 



where : 

1 = ship length 

u^ * operating point speed (15 Kts*1.639 ft/sec/Kt) 



Computer Program #2 is the DSL program which models 
this system. The curves of figure II-5 and II-6 exhibit the 
responses of various step rudder commands. These are 
tabulated and cross referenced in table II-4. 

These responses show the characteristics of a realistic 
rudder in that the rudder is never allowed tc slam into the 
steps. lhey exhibit the time delay between command and 
response which is a function of the rate of response (2.0 
degrees/sec) . A control system design with this scheme is a 
much more difficult problem than one with an idealized 
rudder (step response) because the entire rudder control 
system beccmes quite nonlinear. 
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TABLE II-4 



Piqure 

IZ-5 
II-5 
II- 5 
II-5 
II-5 
II- 5 
II-6 
II-6 
II-6 
II-6 
II-6 
II-6 





BODDER COMMAND AND 


RESPONSE 


Curve 


Rudder 
Command (derj) 


Initial 


1 


+ 30.0 


-30.0 


2 


+ 25.0 


-25.0 


3 


+ 20.0 


-20.0 


4 


+ 15.0 


-15.0 


5 


+ 10.0 


-10.0 


6 


+ 5.0 


- 5.0 


1 


+ 30. 0 


0.0 


2 


+ 25.0 


0.0 


3 


+ 20.0 


0.0 


4 


+ 15.0 


0.0 


5 


+ 10.0 


0.0 


6 


+ 5.0 


0.0 



A 



I 




A 
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Ti<^e (sec) 



Figure II-6 
Rudder Responses 
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C. ENGINE RESPONSE 



Figure II-7 portrays a complex model of a gas turbine 
propulsion plant[ 5 ]£ 6 ]- This model contains the elements 
required for a complete dynamic study of the system. For 
the purpose of this thesis, such a complicated model is not 
required if the overall input-output relationship can be 
established. 

Reference 5 establishes an • output speed (0) 
relationship for a step input of desired speed (Ud) and is 
redrawn as figure II— 8. The relationship appears to be that 
of a first or second order system with a time delay. 



The system equations for a first 
with a time delay nay be written as: 



spun 

SEBEBE 

SBEOOT 



-Ts 

= K*SPDDES*e 
= (SEDIN-SPDCUT) *G 
5PDEBB dt 



-f 



Which yeilds the transfer function: 



-Ts 

SEECOT(S) G*K*e 



SPEIN(S) S+G 



order approximation 



Which is block diagrammed in figure II-9. 







Figure II-9 

Propulsion Plant Low Order Model Block Diagram 
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Figure II-7 

Turbine Propulsion Plant Block Diagram 
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Propulsion Plant High Order Model Step Response 



From figure II-8, the time delay, system gain, and time 
constant can be estimated as: 

T = 4.88 sec 
K = 0.9877 
G = 0.092 

Computer Program #3 was used to obtain the step input 
response. The original complex system output and the low 
order approximation are compared in figure 11-10. As 
indicated in this figure, the two responses are very close. 
Considering the linearized approximations made in the 
eguations cf motion, this response is accurate enough for 
system study use and is used as the model fcr speed control 
in chapter III. 

Similar methods may be used to obtain simplified 
low-order models for other- high-order propulsion systems now 
in use (e.g. pressure fired boiler systems, 1200 lb. 
systems, etc.). They may not, however, simplify to a first 
order approximation suitable for system study. A method of 
computer determination of low-order models of high-crder 
systems is contained in ref. 7 and may be mechanized for 
this purpcse. 
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Figure 11-10 

Stsp Response Comparison of Low and High Order Propulsion 

, Plant Models 



D. EXTEBM1 FOBCES 



The modeling of ship dynamics cannot be complete 

without the introduction of external forces which perturb 

its responses. These forces are caused by many factors and 

some are mere relevant than others in the scope of this 

thesis. The two that are considered can cause substantial 
♦ 

perturbations that must be modeled and eventually accounted 
for in the ccntrcl system design. 

1 . Two Sh ip s in Pr oximity 



Whenever two ships operate in close proximity (less 
than 250 feet) , suction and pressure forces between hulls 
are present. Studies have been conducted on the Mariner 
hull£ 1 ] which have produced data for construction of a 
family of curves fer two ships passing on the same heading. 
No data has teen found for the cases of two ships not cn the 
same heading. Other restrictions on the work presented in 
ref. 1 are that the ships are of the same type and of 
similar hull ratios. 

Interactive effects vary as the sguare of speed. 
However, this is only true if both ships are at the same 
speeds. The interaction modification factor is based cn the 
normalized speed of 15 kts. This factor can thus be written 
as: 



2 

SPDE = CDOT1 

Exact effects on the interaction forces and moments in the 
situation where the ship's speeds differ are not available. 
This is inconvenient since during the approach phase, the 
normalized speed of the approach ship (ship B) can be as 
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high as 1.5. If the effect on ship B is as stated above but 
with its own speed causing the interaction modification, the 
interactive forces and moments can be 2.25 times greater 
than without speed considerations. 

Without the ability to pin down this relationship, it 
was ignored in the development of the control laws presented 
in this thesis. Appendix C was written with the expressed 
intent of illustrating the effect of modifying the 
interactive forces and moments to the extremum case 
mentioned above. It must be realized that this case is not 
considered likely in that it is felt that the interactive 
forces and moments modification on ship B are more apt to be 
caused by the speed of ship A. If this is so, since ship A 
is kept at a constant 15 kts., the interactions need no 
modification for speed consideration in this thesis. 

Beference 1 also gives a method of modifying the 
interactive forces and moments based on different ship 
lengths. For ease of computation and graphical 
presentation, the two ships were considered of equal 
lengths. Tc modify this to ships of dissimilar lengths, the 
resulting bydrcdynamic derivatives must be modified as shown 
in ref. 1 (also shown in appendix C) . 

Since no closed form expression existed for these 
forces, the family of curves reproduced as figures 11-11 and 
11-12 were quantized in subroutine SLOPES (an adaptation of 
the subroutine of the same name from ref. 11). (Appendix C 
contains a curve fitted subroutine that was compiled after 
the completion of the research on this thesis. It was not 
used for any design cr simulation runs except for those 
presented in that same appendix.) An interpolation 
algorithm is used to approximate the intermediate values 
between quantized values and between the curves of the 
family . 
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* 

Note : 



Lateral Distance DY (Ft) 

I 1 I I I 1 

-400 -200 0 200 400 

To convert to normalized lateral distance DY - 
Divide by the ship length L. 



Figure 11-11 

Family of Interactive Y force Curves 



45 




-400 -200 0 200 . 400 

jc 

Note: To convert to normalized lateral distance DY - 

Divide by the ship length L. 



Figure 11-12 

Family of Interactive N Moment Curves 
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The main purpose of this subroutine is to compute the 
interactive forces between ships in the replenishment at sea 
situation anc output the values for perturbation of the 
control ship only. Even though both ships are affected by 
these perturbations, a one ship control system which is 
effective regardless of the other ships motion is 
considered. Consequently, the interactive forces on the 
second ship are ignored. 

Subroutine SLOPES is contained in appendix k. figure 
11-13 is a geographic plot of the two ships passing at 105.6 
feet with their rudders amidships (0 degrees). Figure 11-14 
and figure 11-15 show the magnitude of the lateral force (Y 
force) and rotational moment (N moment) of the reference 
ship on the ship making its approach (control ship) . The 
reference ship is at 15 kts. and the approach ship is at 
22.5 kts. The control ship starts its approach 5 ship 
lengths (2639.0 feetl astern and 0.4 ship lengths (211.12 
feet) laterally displaced. The most graphic portrait cf the 
effects cf these forces and moments appears in the yaw 
changes which are presented in figure 11-16. From these 
figures it becomes readily apparent that these perturbations 
cause viclert motions of the ship which must be accounted 
for in any control system development. Throughout the 
development of such a control system in chapter III, these 
forces and moments are considered inherent in the model for 
HAS control. 
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_2. Waves 

The modeling of sea state in the form of waves and wave 
interactions has occupied the time of many researchers 
[•][*]. The exact formulation of waves will net be 
accomplished in this thesis. Since the main concern here 
will be to test the control scheme developed in chapter III, 
a much simplified wave generator can be used. To introduce 
the reguired experimental perturbations on the designed 
control system a periodic wave system with a fundamental 
freguency ard its second harmonic is used. Some small 
random wave properties are introduced that ride on these two 
sinusoids. A simple expression of this combined wave can be 
written as: 

2 

S = WF*sin(WE) + (PI»WF /WL) •sin(2 # WE) +WF«WRV»sin{WE) 
where : 

W denotes the Wave 

WF is the Wave Force 

WE is the Wave Encounter radian Frequency 

WL is the Wave Length 

WHV is the Wave Random Variable 

PI is *. 1415926 

With this wave as a basis, a method of modeling this in 
the dynamic environment of the total ship simulation was 
defined. The modeling includes the introduction of this 
wave into the three degree of freedom equations of motion. 
To accomplish this a set of defining relationships were 
developed. First the general wave direction is input to 
establish the direction of the wave encounter on the ship. 
If the ship direction is IAWDP2 and the wave direction is 
RD, the expected wave direction is defined as: 

EWC = WD-IAWDP2 
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Next the wave encounter frequency (radian frequency) can be 
established with knowledge of the ships normalized true 
speed (CDCT2) , wave length (WL) and normalized wave velocity 
(WV) . The wave encounter frequency (WEF) is then: 

WEF = 2«PJ> (CD0T2+WV* cos (EHD) )/WL 

The total wave encounter (WE) is nothing more than the 
wave encounter frequency (WEF) times time. This gives the 
wave encounter radian frequency required in the simple 
expression for the combined wave previously shown. 

This dees not complete the task, since the individual 
wave forces of each degree of freedom must be derived for 

this general wave expression, namely the components of WF. 

Again a much simplified version of the more complex real 
life wave forces were used. The X and Y forces are 

considered first. These can be modeled as cosine and sine 
functions of the expected wave direction (EWE) such that: 

WFX = WF*cos (EWD) 

WIS = WF«sin (EWD) 

where WF is the total wave force of the encountered wave. 

The rotational N forces are a little more difficult to 
establish. Ey considering that no rotational forces are 
created by a wave directly on the bow or stern, or directly 

off the beam, and that it is maximum when the wave is at 45 

degrees off the bow or stern, a much simplified 

approximation is developed. Realizing that this method is 
very crude, the N force can be written as: 

WF N = WF«sin (2«EWD) 

To add mere creditability to the wave mcdel, a random 
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variable is added to the wave force at the waves fundamental 
frequency. 2 gaussian (normal) distribution was chosen with 
a zero mean and a standard deviation of one-tenth the 
maximum allowable force of WF. A zero mean signifies that 
the expected amplitude of the random wave is 0.0, while the 
standard deviation signifies that 68% of the random waves 
will have amplitudes less than one-tenth of the maximum 
allowable force of WF. Also, 94% will have amplitudes less 
than one-half the maximum allowable force of WF. This small 
added perturbation allows for verification of the model 
simulation with a stochastic force, which in turn adds 
creditability to the developed ccntrol systems. 

What remains is to define the total wave force (WF) . 
It is important not to fall into a common simulation pitfall 
which inevitably causes unneeded design changes. A sea 
state does net increase at an infinite rate. It therefore 
is incorrect to start a simulation with initial conditions 
set for calm sea and immediately introduce a high sea state 
perturbation. The initial large perturbation transient can 
give results that are not only unrealistic, but can cause 
the model and control system to produ3e unstable results. 
This is especially true in this case since the linearized 
(small perturbations about an operating point) equations of 
motion are used. 

With this in mind, a ramp feed in of the wave force 
with a limiter at the desired maximum wave force ( WFM A) was 
used. The slope of the ramp was established to impart 
minimum initial transients, yet increase the wave force to 
an acceptable testing level within the time frames of the 
simulations cf chapter III. The slope is designed such that 
the maximum wave force is reached in 94.815 seconds actual 
simulation time (4.548 seconds problem time). 

Computer program #4 was used to verify the wave action 



50 






I 






model. Table II-5 on page 56 indicates the figures produced 
and changes in input wave length (WL) and wave direction 
(WD) for each run. The input parameters that were constant 
for all runs are tabulated below: 

YAHEE2 = 0,. 0 
CECT2 =1.5 
* 

8S = 5.0 
WFBA = 0.1137 

* 

NCTE ; WS is the unnormalized wave speed. Conversion to 

normalized wave velocity is: 

HU = WS/15.0 

Introduction of the wave forces is accomplished by 
multiplying the established wave forces by the rudder 
hydrodynamic coefficients for the individual reference 
directions. This effectively scales the wave forces to the 
ship model being used. The wave force result is coded in 
the ship simulation program as follows: 

IF1 2 = FA1*D2+YY2+KA1*WY 
IF 22 = KE1«D2+YN2+KB1*HN 
IF32 = KC1«D2+NC2+KC1«HX 

Detailed results of the wave force effects are given in 
chapter III and will not be dealt with here. 
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Interactive Forces Effect on the Geographic Plot 
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Figure 11-14 
Lateral (Y) For 
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Figure 11-15 
Rotational (N) Moments 
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Figure 11-16 

Interactive Forces Effect on Yaw of the Control Ship 



Run 


* 

Figure 


Input Pc 
* 

WL 


irameters 

* 

WD 


1 


11-17 


0.5 


015 


2 


11-18 


1.0 


015 


3 


11-19 


1.5 


015 


4 


11-20 


0.5 


030 


5 


11-21 


1.0 


030 


6 


11-22 


1.5 


030 



* 

NOTE: WL is given in ship lengths 

WD is given in degrees 

Curve numbers of all runs corresponding to 
wave force components are : 

Force Curve 

WX 1 

WY 2 

WN 3 



Table II-5 

Wave Simulation Listing 
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t Figure 11-21 
Wave Simulation Run #5 
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III. REPLENISHMENT AT SEA 



A. HEADING CONTROL 
1 . Contro l Ch oice 

Many studies involving replenishment at sea (RAS) have 
treated the problem as a multivariable systemf 1 0 ][ 1 1 ][ 12 ] 
[ 13 ]. Academically, there is nothing wrcng with this 
approach. However, as a practical system it leaves much to 
be desired. The key drawback in the multivariable system is 
the inescapable dependency on a command and control link 
between the replenishment ships. The unreliability of UHF 
communications at these close distances is a much 
experienced phenomenon!. It is felt that any knowledgeable 
commanding officer would not entrust the safety of his ship 
to such a guestionable link. An alternate method which is 
described here is a modern extension of the long trusted 
"seaman's eje" concept, where the sensors and control 
devices must be self contained on the ship making the 
approach (hereafter referred to as the receiving ship or 
ship B) . 

In all present day RAS operations, the ship on which 
the approach is being made (hereafter referred to as the 
supplying ship or ship A), must maintain the replenishment 
course and speed. The receiving ship accomplishes the 
maneuvers tc maintain station relative to the supplying 
ship. 



The parameters which are presently measured "visually" 
are relative position (in both the X and the Y directions) , 
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relative heac (usually in reference to ordered replenishment 
course) , and relative motion in the X direction (for speed 
matching) . These parameters are usually visualized by the 
ccnning officer who in turn gives corrective orders to the 
helmsman. The helmsman must then translate these verbal 
orders into rudder and speed commands through the helm and 
lee helm ccnscles. The accuracy of the execution of the 
conning cfficer*s orders is extremely dependent cn the 
ability cf the helmsman and throttleman. This system can be 
guite effective, and it can also be guite disastrous. This 
fast reacting and constantly changing environment lends 
itself tc breakdown in communications and manifests the 
inability of seme individuals to cope with the reguired 
critical man-machine interfaces. 

To eliminate these problems, present state of the art 
digital computers and sensors are available for immediate 
implementation of a completely automatic ship control 
system. Such a control system may be installed on 
individual ships and be used for RAS without the reguirement 
of having the matching installation on the ether ship cf the 
hcckup (another drawback of the multivariable approach) . 

2_. Ccnt rcl Method 



One cf the many pitfalls that may be encountered in 
digital simulation is the reality of the parameters that are 
measurable iE the real world situation as opposed to those 
that are incidently available in the simulation. With this 
fact as a keynote. Subroutine RBMEAS (Range and Bearing 
MEASurement) was developed. This subroutine, as listed in 
appendix l, defines the forward (FWD) and after (AFT) 
relative and true bearings, and ranges from the receiving 
ship to similar points on the supplying ship. Figure III-1 
delineates the terms used in the subroutine. The SDFn terms 
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Figure III-l 
Measurement Techniques 
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are used tc position the bow and stern sensors and 
reflectors cn ship B and A, respectively, in geographic 
coordinates as a function of ships head. The FWD distance 
on the X coordinate is ADFX and the Y coordinate is ADFY . 
Similarly, the AFT distances are ADAX and AEAY. R1 and R2 
are the FWE and AFT ranges measured by a highly accurate 
ranging device installed on ship B. This same ranging 
device, if properly provided with a pinpoint reflector on 
the supply ship (ship A), will give accurate relative 
bearings FWD and AFT., B1 and B2 respectively. The distance 
between sensors may be varied, but as a rule should be kept 
as far apart as possible to allow maximum sensitivity. The 
distance used in this thesis is 1.0 (one ship length), and 
the distances were considered the same for both ships. This 
is net a necessary condition and may be changed to suit the 
situation . 

Subroutine RBKEAS assumes highly accurate sensors in 
both ranee and tearing measuring ability. Such sensors are 
presently available in the form of Radar alt imeters[ 1 ♦ ] and 
Laser ranging devices. Another possibility for a measuring 
method is a single sensor time sharing to obtain range and 
bearing tc both reflectors from a single device. Such a 
single sensor scheme is sketched in figure III-2. 

Once the FWD and AFT parameters are available, they may 
then be used to determine other desired quantities. 
Subroutine HDGRAS (HeaDinG control for RAS) was developed to 
output the desired heading corrected for heading difference 
of ship A and B and the projected correction for distance 
error. This subroutine is listed in appendix A. The center 
range and hearing are the average of the FWD and AFT range 
and hearing output from Subroutine RBMEAS . 

The additional heading due to distance is projected as 
if ship E maintained its present course until it was 
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perpendicular to the center of ship A. The reasoning behind 
this is illustrated in figure III-3. If the present course 
will cause ship B to arrive on the station desired (DS) , no 
heading chance is required. The expression for PSIADC 
Additiocal heading due to Distance Correction) is: 

PSIADC = ESENS* (DDC+DA-SIN (AA1) ) 

where : 



RSENS = Bange SENSitivity gain 

DDC = Distance Desired Corrected for side of approach 

DA = center Distance Absolute (range) 

AA 1 = 360 degrees - relative bearing of center position 

The beading difference of ship A and B is desired 
since, even if the range when alongside is correct, a large 
disparity in heading cannot be tolerated. It is realized 
that some beading difference (crabbing) is necessary to 
maintain the distance. This crabbing is due entirely to the 
pressure fcrces modeled in chapter II. This heading 
difference is found by computing the difference in the 
perpendicular projection between the PHD and AFT 
measurements and finding the arcsin of this difference 
divided by the distance between sensors. Figure III-4 
indicates a sample of this procedure. 

The expression for total desired heading is given as 
fellows : 

PSIDES = PSIAEC+WTSENS*PSIDIF+PSIE 
where : 

PSIEIF = ? additional heading due to heading DIFference 

WTSENS = BeighTed heading difference SENSitivity gain 
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Figure III-3 
Distance Logic 
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Figure III-4 

Heading Difference Calculation 



70 



PSIDES = $ (heading) DESired 



Throughout the subroutines and main DSL programs, the 
Function EEGBAD (conversion of degrees to radians and 
radians tc degrees) is used freely. An explanation and 
listing cf this function are presented in appendix A. 

The angular velocity of the receiving ship’s head is 
also cf concern in the HAS situation. This quantity may be 
thought of as similar to tachometer feedback in a simple 
serve ccntrcl system; and is necessary to damp out the 
response (the responses of this control system without this 
feedback is presented in the latter section of this 
chapter) . 

The desired rudder command is a combination cf the 
desired heading, angular velocity feedback, and a rudder 
gain as fellows: 

Desired Eudder = (YAWD2-PSIDEE+BD0TFB) *BGN 
where : 

YAHD2 = heading of ship B (in degrees) 

PSIEED =E SIDES (in degrees) 

BDOTEB = VFBG-BECT2D 

VFG = Velocity FeedBack Gain 

BEOT2E = angular velocity of ship B heading angle (in 
degrees/sec) 

EDOTFB = angular velocity FeedBack 

BGN = Budder GaiN 

The convention for rudder response dictates negative 
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rudder as being right rudder, which causes positive yaw. 
This necessitates making the desired rudder the negative of 
the forcing function and feedback quantities. The block 
diagram cf figure III-5 presents the control loop from 
measurement inputs to desired rudder command. 
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3 . Opti mizat ion 



Thus far the control choice has identified four gains 
(BSENS, KISERS, 3GN, VFBG) that must be set for proper 
position attainment. The nonlinear nature of the system 
which appears in the form of distance measurements, 
interactive forces and rudder modeling do not allow for 
straight forward determination of these gains with normal 
optimal control theory. 

a. Techn iq ue 



Grossly nonlinear systems require special handling to 
determine proper gain settings. The method chosen for this 
purpose was an optimization algorithm developed by H. «3. Box 
(programmed locally as subroutine EOXPLX) . This subroutine, 
listed in appendix A, was used to locate the cost function 
saddle point in four dimensional space (the dimensions being 
the previously mentioned gains) . The drawback associated 
with this method is the necessity of iterating the complete 
nonlinear simulation within function FE for every evaluation 
of the cost function. The gains sought were found, but 
unfortunately only after 2 1/2 to 3 hours of CPU time with 
every 400 iterations allowed. 

The mechanics involved in optimizing the chosen cost 
function include reguired sub-calculations in many functions 
and subroutines. Figure III-6 is a flow chart which 
demonstrates the steps, subroutines and functions required. 

Initial optimization was accomplished for one set of 
initial conditions. By looking at the RAS situation, a 
probable set of circumstances were envisioned. The scenario 
setting is the approach phase where the replenishing ship 
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Figure III-6 
Optimization Flow Chart 
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starts a wide approach at 0.4 ship lengths (211 feet) 
lateral displacement and 5 ship lengths (2639.0 feet) astern 
of the supply ship. The desired final position is alongside 
at 0.2. ship lengths (105.6 feet) lateral separation. Eoth 
ships have the same initial heading (YAW angle). The supply 
ship is at 15 kts. (1.0 normalized speed) and the receiving 
ship makes its approach at 22.5 kts (speed control will be 
covered later in this chapter) . 

b . Cost fun ction 

Normal costing of displacement error with the integral 
of the sguared error (ISE) was considered as the 
optimization tool in subroutine BOXPLX. However, this type 
cf performance measure would weigh the initial displacement 
error egually with the final position error. This problem 
can be circumvented by comparing the displacement error to a 
pre-ccmputed reference track instead of to the desired 
displacement. For the envisioned scenario, it was conceived 
that the cost function should weight the distance 
displacement heavier when the ships are alongside than when 
the approach is started 5 ship lengths astern. 

This was accomplished by using the integral cf time 
times the absolute error (ITAE) as the optimization 
performance measure. The reference displacement was 
considered the desired position displacement. The object 
furcticn can then be written as: 




to 



dt 



where 



DD = Desired Distance 

AEY = Actual Displacement in the Y direction 
t = time 
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A performance measure that is designed to obtain good 
performance must also take into account other factors 
besides just position accuracy. Consequently, another cost 
criterion was decided upon which would also set the gains to 
reduce tie amount cf rudder control required when alongside. 
This particular feature is derived from the desire net to 
over control with the rudder in such close proximity to 
another vessel. The inclusion of this term in the 
performance measure is weighted by unity while the distance 
accuracy is weighted by a factor of 10.0. This will tend to 
allow rudder action if the desired position is not 
maintained. The final approach phase cost function for 
obtaining optimum gains has the form: 



CEJ 




1 0 . 0 • | DD- ADY | + 1 . 0* | D2 ) ) 



dt 



where the additional term is: 

D2 = rudder response of the replenishing ship 



c . R esu lts 



In the process of deciding on the best gain definitions 
previously mentioned, many optimization runs were made. 
Each set cf gains were then simulated in a corresponding DSL 
program to obtain performance confirmation. Many of these 
runs did not live up to expectations; causing re-evaluation 
of the control scheme until the one presented in this thesis 
was formulated. 

Table III-1 shows the input upper and lower limits of 
search (ED, EL), starting value guess (XS) , optimum gain 
settings (Output) and associated object function value (CBJ) 
for 20.0 seccnd normalized time simulation run in function 
EE. These values were then introduced into the DSL 
simulation program listed as program #5. The results of 
this simulation are shown in figures III-7 thru III-12. The 
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Approach Phase Optimization Results 
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Figure III-7 

Approach Phase Yaw Result 
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Figure III-8 
Approach Phase Y Forces 
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Figure III-9 

Approach Phase N Moments 
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Figure III-10 

Approach Phase Geographic Plot 
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Figure III-ll 

Approach Phase Lateral Distance DY 
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Figure III-12 

Approach Phase Rudder Response 



indicates 



geographic plot cf figure III-10 indicates excellent 
positioning in the lateral direction while the rudder 
response cf figure III-12 shows that it settles out to a 
fairly constant steady state value as the ship settles into 
its desired position. The time coordinates in all plots are 
shown in actual full scale time. 



d. Control Test ing 



Now that the "proper 1 ' gain settings were obtained, more 
extensive testing of the control system was required. Three 
different tests were contemplated: (1) allow a large 
perturbation turn of the reference ship (supply ship) , (2) 
start approach of the receiving ship (control ship) from 
different initial conditions cf lateral and horizontal 
displacements, and (3) induce external perturbations in the 
form of wave forces. 

The first test was simulated by turning the reference 
ship by normal rudder action of figure III-13. This turn 
with 5 degrees rudder accounted for a total reference yaw 
change of 15 degrees. The rudder action of the controlled 
ship shewn in figure III-14 was as expected. However, the 
distance maictainment portrayed in figure III-15 was totally 
unacceptable. The maximum excursion from the desired 
distance cf 105.56 feet (0.2 normalizes distance) was 55.419 
feet (C.105 normalized distance). Variances of this 
magnitude cannot be tolerated in the RAS environment. 

Faced with this situation, the tact chosen was to 
re-evaluate the gains for the new scenario which is called 
the turn phase. In this phase the initial conditions assume 
steady state positioning alongside such that the lateral 
position displacement (DY) is equal to the desired distance 
[105.56 feet (0.2 normalized)] and that the horizontal 
position displacement (DX) is 0.0 (alongside). 
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Turn Phase Rudder Action of Reference Ship 
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Turn Phase Rudder Response 
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Figure III-15 

•Turn Phase Lateral Distance DY 



Seme initial perturbation is introduced by assuming the 
relative yaw angle when alongside is negligible. 

e_. second optimiz ation 

The same pxocedure was followed in obtaining gains that 
would optimize a chosen cost function. Figure III-6 still 
applies except that function FE A is replaced by functicn FEB 
(listed in appendix A) to simulate the new conditions. 

Cost function criteria change in this instance since 
the ships start at the desired position and optimally stay 
at the same relative positions. Also, the rudder response 
to such a large turning perturbation must be free to cause 
achievement of the desired position. Due to these 
considerations, the integral of the absolute error (IAE) 
performance measure was chosen for the optimization 
criterion and can te written as: 



Table III-2 shows the results of the turn phase 
optimization and the comparison with the approach phase 
gains. Again CSL simulation was performed using the turn 
phase scenario. Figures III-16 thru III-21 portray the 
graphical results. The rudder response cf figure III-21 
indicates very sensitive response to the interactive forces 
shewn in figures III-17 and III-18. The lateral distance 
separation cf figure III-20 indicates excellent position 
maintainment with maximum excursion error of only 2 feet 
(0.0038 normalized). This minimal variation is well within 
that which can be tolerated in the RAS environment. 




dt 
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Figure III-16 
Turn -Phase Yaw Response 
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Figure III-17 
Turn Phase Y Forces 
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Turn Phase 
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Figure III-19 

Turn Phase Geographic Plot 
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Phase Lateral Distance DY 



Roooe« Response 
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Turn Phase Rudder Response 



f. Ccntirue d Cont rol Testin g 



To alleviate suspicions that the response from the 
gains obtained in the approach phase could be improved by 
the gains obtained in the turn phase, a simulation of the 
approach phase was accomplished with the new gains. Figure 
III-22 is the graphic display of the effect of these gains 
on the approach phase lateral distance positioning. 

i 

Careful analysis of the results thus far clearly 
indicate the need for an adaptive control scheme to allow 
gain adaptation to meet the design specifications. A full 
adaptive control scheme for systems of this type is outside 
the scope of this thesis. References 15 thru 22 are 
indicatiors of some of the literature available for pursuit 
of a completely adaptive control system. 

What was done here is develcp-ment of a simple algorithm 
to sense when the conditions were adequate to switch from 
one set cf gains to another. This may be done with the two 
sets of gains developed thus far. However, for the sake of 
simulation efficiency, a third set of gains was introduced. 
This third set amounts to a change of one approach gain 
(RSENS) which has previously been defined as the -range 
sensitivity gain. The simulation efficiency is increased by 
decreasing the time required for the approach phase to reach 
steady state. A consequence of this procedure is a 
reenforcement of the need for a completely adaptive control 
scheme. 

Repeated simulation revealed that commencing the turn 
(in effect switching gains) , before a reasonable steady 
state was reached caused results similar to those shewn in 
figure III-22. An increase in RSENS to a value of 4.0 when 
the lateral separation error is less than 0.05 (normalized) 
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Figure III-22 

Approach Phase Lateral Distance DY 



and greater than 0.005 (normalized) forces acceptable steady 
state in approximately 1/2 the time previously required 
osirg a single set of .approach gains. 

Subroutine SJiTCH (listed in appendix A) incorporates 
this simple adaptive gain schedule with a counter mechanism 
to sense when steady state is reached. Further study 
indicated a need to damp the yaw oscillations to a g-reater 
extent if the y^w velocity ( BD0T2D) exceeded 2.0 degrees/sec 
when the gaits are initially switched to the more sensitive 
ones of the turn phase. This is an artifical adaptive gain 
for VFEG caused by computer time restrictions prevalent in a 
full scale computer simulation where both the approach and 
turn phases are desired. If the gain switching point is 
moved up in time, as would normally be the case in a real 
life situation, this damping increase would not be required. 

The results of the full scale simulation using computer 
program #6 are shown in figures III-23 thru III-34. The 
approach phase plots, of figures III-23 thru III-28 show 
definite improvement over that previously shown in figures 
III-7 thru III- 12. Figure III-27 indicates that the 
overshoot is reduced to 10.6 feet (0.02 normalized) as 
opposed to 17.9 feet (0.034 normalized) that was prevalent 
in figure ill- 11, 

The turn phase plots are shown in figures III-29 thru 
III-34 and shew responses very similar to those shewn 
previously in figures III-16 thru III-21. The only 
significant differences occur in the initial responses which 
are due to the incorrect initialization when the turn phase 
was simulated individually. 
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Figure III-23 

Approach Phase Yaw Response 
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Approach Phase 
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Figure III-25 
Approach Phase N Moments 
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Figure III-26 

Approach Phase Geographic Plot 
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Figure III-27 

Approach Phase Lateral Distance DY 
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Figure III-28 

Approach Phase Rudder Response 
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Figure III-29 
Turn Phase Yaw Response 
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Figure III-30 
Turn Phase Y Forces 
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Turn Phase 
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Figure III-32 

Turn Phase Geographic Plot 
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Figure III-33 

Turn Phase Lateral Distance DY 
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_g_. Varyin g Initia l Condition s 



The results obtained in the previous section are most 
gratifying tut actually incomplete. This system must work 
for ether initial- conditions quite different from those 
envisioned in the optimization scenario. The initial 
approach can realistically commence at points other than 5 
ship lengths astern and displaced by 0.4 ship lengths. 

Ey simulating this system with varying initial 
positions, the relative efficiency and worth of the control 
system can he observed. This was done in successive test 
runs whose initial conditions and corresponding plot figures 
are tabulated in table III-3. For the sake cf brevity only 
these figures required to illustrate the relative efficiency 
of the ccnticl system are included. The corresponding 
initial optimization simulation figures are listed for cross 
reference. The turn phase plots for all runs except 4 and 6 
exactly match that of the initial simulation and are not 
repeated here. 

Buns 3, 5 and 6 were accomplished to show that no 
ambiguities exist in the control scheme to prohibit adeguate 
real life initial conditions. Run 3 simulates the situation 
most often encountered by this author in the RAS 
environment. This scenario starts the control ship dead 
astern at 5.0 ship lengths and brings ix alongside at 0.2 
ship lengths lateral separation. 

Run 5 is a situation where the approaching ship is 
purposely placed out of position on the wrong side for 
approach. The control scheme adequately corrects the 
placement error and will do so for all cases of this type, 
provided that there is adequate maneuvering room astern of 
the reference vessel (in this case 2.6 ship lengths was 
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Initial Condition Simulation Cross Reference 
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Figure III-35 

Approach Phase Run #1 Yaw Response 
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Figure III-36 

Approach Phase Run #1 Geographic Plot 




116 



Figure III-37 

Approach Phase Run #1 Rudder Response 
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Figure III-38 

Approach Phase Run #2 Yaw Response 




118 



Figure III-39 

Approach Phase Run #2 Geographic Plot 
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Figure III-40 

Approach Phase Run #2 Rudder Response 
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Figure III-41 

Approach Phase Run #3 Yaw Response 
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Figure III-42 

Approach Phase Run #3 Geographic Plot 
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Figure III-43 

Approach Phase Run #3 Rudder Response 
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Figure III-44 

Approach Phase Run #4 Yaw Response 
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Figure III-45 

Approach Phase Run #4 Geographic Plot 
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Figure III-47 

Approach Phase Run #5 Yaw Response 
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Approach Phase Run #5 Geographic Plot 
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Figure III-49 

Approach Phase Run #5 Rudder Response 
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Approach Phase Run #6 Yaw Response 
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Figure III-51 

Approach Phase Run #6 Geographic Plot 
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Figure III-52 

Approach Phase Run #6 Rudder Response 
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Turn Phase Run #4 Yaw Response 
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Figure III-54 

Turn Phase Run #4 Geographic Plot 
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Turn Phase Run #4 Rudder Response 



X 




135 



Figure III-56 

Turn Phase Run #6 Yaw Response 
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Figure III-57 

Turn Phase Run #6 Geographic Plot 
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Figure III-58 

Turn Phase Run #6 Rudder Response 



experienced which gives 1.6 ship lengths bow to stern 
clearance) . 

The purpose c£ run 6 is to provide simulation for an 
approach from the opposite side again disproving any concern 
for ambiguity in the trignometric measurement scheme 
utilized. In all runs it must be emphasized that DD is the 
positive absolute distance desired and that IS provides the 
code flag for the desired side of approach. The system will 
work with DD set to some negative quantity; but the side of 
approach will reverse itself and the position placement will 
be correct, but on the side not desired. 

Hun 4 takes the desired distance in to 0.15 ship 
lengths (80.0 feet). This distance is usually the minimum 
desired by a prudent seaman. Again, even with this minimum 
distance, the control system performs up to desired 
standards. The importance of this, run cannot be overlooked. 
Performance of the system at this extremum indicates that 
the gains utilized are correct for all expected conditions 
encountered in calm seas. Figures III-59 thru III-64 
portray the remaining plots obtained in run 4. 

h. Perfo rmance in S ea State 



The calm sea performance of the heading control system 
is only part of the system testing reguired. Of even 
greater concern is the adequacy of the control when sea 
state is introduced. Section D.2. of chapter II models the 
three components of waves with two sinusoids and a small 
random impulse wave. These forces were introduced into the 
total RAS simulation as shown in computer program #7. In 
this program the wave length (WL) is set to one ship length 
and the wave direction (HD) is -015 degrees true. This 
scenario allows for a port turn into the prevailing sea as 
is common practice in experienced RAS evolutions. By 
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Figure III-59 

Approach Phase Run #4 Lateral Y Forces 
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Figure III-60 

Approach Phase Run #4 Rotational N Moments 
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Figure III-61 

Approach Phase Run #4 Lateral Distance DY 
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Figure III-62 

Turn Phase Run tt4 Lateral Y Forces 
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Figure III-63 

Turn Phase Run #4 Rotational N Moments 



V 




Q 

<D 

O 

g 

fd 

-P 

•H 

Q 

■— I 

vd fd 
I P 

h a) 
H -P 
h fd 
PI 
CD 

u 

3 ^ 

tjl 4*= 

Cm G 

G 

<D 

cn 

td 

x; 

a. 

G 

U 

G 

Eh 



144 



minimizing the perturbation forces on yaw and lateral 
direction, a smoother RAS can be accomplished thus aiding 
safety and comfort during the actual transfer. The wave 
force maximum is taken as 0.05685. Runs were simulated 
which used maximum wave forces in the range 0.1137 to 
0.05685, wave lengths from 0.5 to 1.5 ship lengths and wave 
directions 015 to -015 degrees off the initial replenishment 
course. The control system handled all of the perturbations 
well except for the cases of a wave length of 1.5. This 
length of wave with a force of 0.05685 exceeded the control 
systems capability in that the steady state conditions were 
not met before a turn was commenced. Figure III-65 shows 
this instability in the lateral distance DY of the turn 
phase. It is felt that the modeling inadequacies of the sea 
state development of chapter II coupled with a simple 
adaptive gain scheme are the source of the problem. This 
sane phenomenon is covered in greater detail in the 
longitudinal position offset testing portion of the velocity 
control section of this chapter. 

Problems of this type also manifest themselves in some 
cases when the wave force maximum (WFM A) was close to the 
0.1137 value. If the sea state becomes excessive, which 
this value represents, a different gain schedule or, at 
best, a mere complex adaptive gain scheme is called for. 

The plots produced by computer program #7 are presented 
as a representative indication of the effectiveness of the 
control system in the presence of a sea state. Figure 
III-66 gives the yaw results of the approach phase which 
indicates the effect of the wave action. The corresponding 
rudder action of figure III-67 compensates to give the 
smooth lateral distance shown in figure III-68. The wave 
profile is shown in figure III-69 with curve 2 being WY and 
curve 3 being WN. Curve 1 is the WX profile which was not 
used in this run but will be utilized in the speed control 
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section later in this chapter. Similar curves are portrayed 
for the turn phase. Figure III-70 is the yaw difference 
between the two ships (remembering that the reference ship 
is net being perturbed by the interaction forces or the wave 
forces) . Figure III-71 is the lateral distance DY 
maintained by the rudder response of figure III-72. The 
maximum lateral separation . in the turn phase is 0.0037 ship 
lengths (1.S5 feet). The wave profile is shown in figure 
III-73 with the same wave force curve sequence as the 
approach phase. 

Ss can be seen from these plots, the control system 
operates very effectively in the presence of a sea state. 
Again, the development of a much more complex adaptive gain 
scheme is reguired to allow exceptionally high sea state. 
It is felt that the control system presented in this thesis 
is adequate for most situations that are encountered in the 
HAS environment. Only the extreme perturbations that chance 
would allow must be accounted for in a more complex adaptive 
gain scheme. 
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Figure III-65 

Wave Effect on Turn Phase Lateral Distance (DY) WL= 
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Figure III-66 

Wave Effect on Approach Phase Yaw WL= 
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Figure III-67 

Approach Phase Rudder Response to Waves WL= 
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Figure III-68 

Wave Effect on Approach Phase Lateral Distance (DY) WL= 
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Figure III-69 

Approach Phase Wave Profile WL= 
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Wave Effect on Turn Phase Yaw WL- 
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Figure III-72 

Turn Phase Rudder Response to Waves WL- 
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Figure III-73 

Turn Phase Wave Profile WL= 



B. VELOCITY COHTBOL 



One advantage derived from using the linearized 
equations of motion is the decoupling of the velocity 
components from the remaining equations of motion. This 
allows separation of the design procedures for lateral 
separation control and velocity control. Section A. of this 
chapter designed the lateral separation control using the 
simple speed control algorithm shown in figure III-74. This 
control output was used directly as the ship's speed (CDCT2) 
in the model simulation where no attempt was made to use the 
engine response developed in chapter II. Pucction SEDCTB of 
appendix A shews the control used. 
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Figure III- 7 4 
Non-optimum Speed Law 



Because of this decoupling assumption, any valid 
approach speed control can be used, if used consistently. 
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for such a design. However, in the BAS environment, 

complete disasscciation is not possible. Recombination 
occurs in the interactive forces and moments which depend 
upon the longitudinal distance as well as the lateral 
distance. Consequently, speed, which is directly 

responsible fcr the longitudinal distance, has a direct 
relation to the lateral distance attainment and 
maintainuent . 

*• * 

The remaining parts of this chapter deal with the 
development cf a viable speed control algorithm and the 
testing of the designed system. 

1 . Type cf Control 



Hhenever two ships maneuver for replenishment at sea 
(BAS) , the prime considerations are the time required for 
approach and the accuracy of position keeping plus 
conservation of fuel. 

The nonlinear control law of figure 111^75 is designed 
to maintain a preselected approach speed for minimum 
approach time. The proper location of the switching point 
increases the complexity of the solution since the time of 
switching from this speed is determined by the dynamics of 
the nonlinear position attainment loop. Once this position 
is reached, the speed controller is switched dawn to a 
linear pcrticn of the control law to allow control for 
perturbations about the operating position. However, small 
perturbations about this operating point can be tolerated 
and, in fact, are desired to allow for conservation of fuel. 
Selection of this dead zone is wholly dependent on the 
accuracy required for final position. Figure III-75 
indicates a dead zone extending to ±0.001 normalized 
distance which in this case translates to ±0.53 feet. 
Systems fcr which fuel considerations are not a motivating 
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factor nay be designed without this part of the control law 
to allow finer tracking in the position loop- 



SPEED 




The speed, control law as explained above is shewn in 
figure III— 75 for an initial approach speed of SPD02 and a 
final estimated reference speed of SPD01, with ADX being the 
dynamic position feedback defined as the longitudinal 
distance between centers of the ships referenced to the 
controlled ship’s heading- Analytically,, the linear portion 
of the control law is written asr 

SPICTH a — ADX*- (S PDC2— 5 ED01 ) +SPD0 1 

Symetric continuation of the control law accounts for 
operation on both sides of the operating point. 
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2. Optimization 



Using this suck simplified model of chapter H and the 
basic control la a of figure III-75, the desired snitching 
curve can be established*. An optimization subroutine such 
as; subroutine BOXPLX can be used, to iterativelj obtain the 
optimum snitching position (Stf )- for representative initial 
approach speeds., Figure- III-76 is a flow chart of the 
subroutines and' functions reguired for speed control 
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Figure IIX-76 
Optimization Flow Chart 



optimization. The major merit of this nonlinear control law 
stems from the predetermination of the switching point for 
all possible conditions of initial speeds. This apriori 
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knowlege allows for offline computation of the switching 
position prior to commencing the approach. The cost 
function used for optimization is the ITAE which 
accomplishes two ob-jectives. First, it forces the approach 
to be accomplished in minimum time. Secondly, it insures 
that the fuel expenditure will be optimized in the 
elimination cf most overshoot and bang-bang control in the 
dead zone portion of the control law. The final value of 
the position error must be within the specified dead zone 
and the terminal speed must match the reference speed 
(SPD01) . The cost function has the following form: 



J 




(t • 1 ADX | ) dt 



Table III-4 is a comparison of the optimization runs 
with various initial speeds. The values shown for SH must 
be multiplied by the speed differential (SPD02-SPE0 1) to 
obtain the corresponding value of ADX. The max/min values 
show the band cf values which produce the optimum cost. 
This range cf values is attributed to the integration step 
size used in the optimization program. Experience with this 
particular optimization program indicates that erroneous 
values of the switching point are found if the step si 2 e is 
net carefully chosen. The step size may be adeguate for 
integration, but net for location of the switching point. 



The points obtained from the optimization runs are 
plotted in figure III-77. These points define the nonlinear 
switching curve which must be stored in the computer to 
insure optimal operation of the speed control for all 
approach speeds. From here there are many procedure options 
open. These options have as a goal some usable form for 
predicting the optimal switching point for any set of 
initial conditions. One may choose linear straight line 
segments with an interpolation routine, or a closed form 
switching curve polynomial. Due to the availability of a 
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II 


VITIAL Cl 


URVE POINTS 


SPD02 


SPD01 


SW 

MAX 


SW 

MIN 


SW 


COST 






★ 


* 




* ★ 


1.1 




.545 


.545 


.545 


22.340515 


1.2 




.58705 


.58424 


.585 


5.733367 


1.3 




.62656 


.6256 


.626 


2.700768 


1.4 




.6845 


.68234 


.683 


1.672599 


1.5 


1.0 


.73169 


.7283 


.729 


1.223071 


1.6 




.7644 


.76142 


.763 


0.992283 


1.7 




.7945 


.7926 


.7936 


0.861552 


H 

• 

00 




.82178 


.81945 


.82 


0.774621 


1.9 




.8501 


. 8439 


.85 


0.757244 


2.0 




.8673 


. 86375 


.865 


0.730168 


CURVE CHECK POIN' 


rs 


1.5 


1.1 


.6859 


.6823 


.683 


1.668055 


1.5 


1.2 


. 6307 


.6297 


.6302 


2.691659 


1.6 


1.2 


.67965 


.67906 


.6793 


1.659956 



★ 

cpu usage over 4 min. - run not complete. 
** 

cost function based on 20 min problem time 
all others based on 10 min problem time. 

Table III-4 
Optimization Results 



161 




162 



Figure III-77 

Switching Curve Minimization Results 



hybrid configured XDS 9300 digital computer and AGT-10 
graphics terninal, the latter course was taken. 

A polynomial curve fitting algorithm was used to obtain 
the required polynomial coefficients of best fit. This was 
done for polynomials of order 1 thru 5. Th-e coefficients 
and the sum cf the sguares of deviation frcm the original 
points are tabulated in table III-5. The selection cf the 
order to te used is highly dependent on the degree of 
accuracy reguired. In the BAS problem, the average error 
introduced for a first order fit is 8.0 feet(1.07 sec), 

while the fifth order fit introduces an average error of 
1.35 feet (C. 180 sec). Prior acceptance of errors introduced 
by an integration (and problem) step size cf 0.8 sec allows 
for use cf a second order fit without any degradation of 
simulation accuracy [second order average error is 2.848 
feet (0.38 sec) ], The graphic display of figure III-78 
indicates very little difference in the switching curves for 
second tc fifth order polynomial fits. For the sake of 
accuracy, and owing to the computer control methods of this 
thesis, the fifth order polynomial fit shown separately in 
figure III-79 is used for determination of the switching 
point location. 

3. Control Testin g 



A true test of the control law is accomplished when it 
is introduced in a computer program fcr a complete RAS 
simulation. Considering the performance of this controller 
in a complex environment of full scale BAS simulation allows 
maximum verification of the controller design. 

The scenario for this simulation initially positions 
the ships such that the ship being controlled starts an 
approach 5 ship lengths (2639 feet) behind the reference 
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Figure III-78 

First Thru Fifth Order Polynomial Curve Fit Results 
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Figure .111-79 

Fifth Order Polynomial Curve Fit 



ship and displaced 0.4 ship lengths (211 feet) to the right. 
The desired final position is alongside and displaced 0.2 
ship lengths (106 feet) . The heading control system used is 
developed in section A of this chapter. 

The approach phase is accomplished with the speed 
desired and speed acguired shown in figure III-80 with the 
corresponding position attainment exhibited in figure 
III-81 . These plots show excellent switching and optimal 
position attainment. 

The next step is to insure that the position keeping 
loop will maintain the desired position with an induced 
perturbation. This is accomplished by turning the reference 
ship away frcm the control ship a total of 15 degrees to 
observe the reaction of the speed control loop. The 
reference ship's turn causes the relative motion between the 
ships to be altered, making the control ship lag the desired 
position. The nonlinear control system is designed to 
correct this situation as soon as the actual position is 
outside the limits of the dead zone. Figure III-82 displays 
the desired speed and acguired speed for the control ship. 
Figure III-83 indicates that the corresponding position 
deviates frcm the desired by 0.0154 ship lengths (8.13 feet) 
at the maximum excursion. This is well within £he limits of 
acceptability for such a drastic perturbation. 

The introduction of velocity control was accomplished 
by combining the simplified engine response of chapter II 
and the speed control law developed here. By setting the 
speed desired (SPDDES) eguai to the output of Function 
SPC EEC and scaling the speed error (SPEERR) to the 
nondimensional eguations of motion, the velocity loop is 
initiated. The auxiliary eguations added to those presented 
in chapter II are: 



167 






Figure III-80 

v 

RAS Speed Control Approach Phase 
Speed Desired (1) and Speed Acquired (2) vs Real Time 




168 



) 




169 



i0 -00 20-00 30.00 30.00 tO-OC 60-00 



Figure III-82 

RAS Speed Control Turn Phase 
Speed Desired (1) and Speed Acquired (2) vs Real Time 
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Figure III-83 

RAS Speed Control Turn Phase 
Position Attainment vs Real Time 
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SPEDES = SPDREC (ADX , S PD01 , SPD02 , S W) 

CDCT2 = INTGRL (U02 , SPDEHR*L0C) 

These equations are introduced in computer program #8 to 
produce figures III-80 thru III-83. 

Further system study indicates that the reference ship 
speed must be known to a fairly high degree of accuracy. 
Without apriori kncwlege of the reference ship speed, a 
constant bias is introduced. The amount of bias allowable 
defines the permissible uncertainty in the reference ship's 
initial speed. This bias can amount to as much as 0. 1 ship 
lengths (84.48 feet) for a reference speed inaccuracy cf 2.5 
knots (0.1 normalised speed) . However, it is felt that the 
reference speed in any practical situation will initially be 
known to within 0.5 knots (0.02 normalized speed). This 
mere practical error will introduce a bias of only 16 feet. 

Other feedback parameters can be used to offset the 
lack of apriori kncwlege of the reference ship speed. Since 
the reference ship is tracked with a high accuracy range and 
bearing device and the controlled ship's speed is measured, 
a decoupled multivariable scheme is used to further refine 
the reference ship speed. With high resolution devices 
presently available^ 14 ], it is estimated that this can be 
dene practically to within 0.05 knots (0.002 normalized 
speed). This would bring the offset bias to 1.6 feet; well 
within previously defined errors introduced by integration 
step size. 
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ii Long itud ina l Po sition O ffset 



Throughout the development of the heading control and 
speed control, the scenario has followed the condition that 
the final position would be longitudinally alongside. 
Although this is a good assumption for ships of the same 
type, it dees not account for BAS station differences for 
different ship types. To alleviate this disparity, function 
SPEBEC was redesigned to allow pre-planned offset condition 
to exist. Function SPDOFC of appendix A is a result of this 
redesign . 

Simulation runs, with a change of the speed control 
function only, resulted in some unstable conditions existing 
in the heading control loop. The cause of this phenomenon 
stems back to the adaptive gain scheme used and the changes 
made to force the control loop to a steady state value prior 
to a turn. Ey using a favorite ploy of experienced ccnning 
officers, this problem is alleviated. The ploy is to take 
the ship alongside and then either drop back to station or 
surge forward to station. This method is accomplished by 
setting the initial offset (XOFS) to 0.0. The final desired 
offset (XCFSE) is stored and not used until the ship is 
settled out alongside. It is subsequently used as shewn in 
the following Fortran code: 

IF (ATIME.GT. 450.0) XOFS = XOFSD 



This method solved the gain transition problem. It did 
not, however, give a completely stable simulation run. 
Onstatle conditions still existed at the end of the turn 
phase. lhis is not suprising, considering the heading 

control optimization method used. The set of gains 

previously found were for the alongside scenario only. 



173 



Different interactive forces and moments at the offset 
position cause these gains to be no longer optimal- 

By relaming the control loop in the heading velocity 
feedback gain (VFBG) , sub-optimal control at all practical 
offset positions is achieved- The gain- VFBG was changed 
from: Q. 064028 to 0-T in the turn, phase adaptive gain 
schedule without significant less of control efficiency for 
alongside operation (2-3 feet maximum excursion vice 2-0 
feet previously obtained) - Subroutine SWTCHF of appendix A- 
reflects the gain change and offset calculations required. 
.Computer program #9 incorporates the changes required for 
offset simulation. Table III-6 is a cross, reference listing 
of the plots obtained. From these figures, the effect of 
different longitudinal positions is readily apparent- An 
offset of 0-1, eguating to 52.3 feet, causes greater lateral 
excursions when astern (XOFSD = -0.1) of the alongside 
position than when ahead (XOFSD = 0.1). The longitudinal 
position maintainment, however, is essentially the same in 
all cases. 
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Note: These plots are the same as those obtained from 

computer program #8 and are not repeated here. 

Table III-6 

Position Offset Testing Cross Reference 



174 



An alternative to the method shown here is again a 
completely adaptive gain scheme which would achieve optimal 
control instead of the sub-optimal control settled for here. 
The alternative may become even more important if the 
nonlinear terms of the equations of motion are considered. 
This would couple the heading and speed control designs to a 
larger estent than encountered in the interactive forces and 
Elements. 
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Approach Phase 
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Approach Phase 
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Approach Phase Run B Lateral Distance 
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Figure III-88 

Approach Phase Run B Speed Response 
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Figure III-91 

Approach Phase Run C Yaw Difference 
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Figure III-92 

Approach Phase Run C Speed Response 
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Figure III-93 

Approach Phase Run C Longitudinal Position DX 
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Figure III-94 

Turn Phase Run A Lateral Distance DY 
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Figure III-95 

Turn Phase Run A Yaw Difference 
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Figure III-96 

Turn Phase Run B Lateral Distance DY 
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Figure III-97 

Turn Phase Run B Yaw Difference 
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Figure III-98 

Turn Phase Run B Longitudinal Position DX 
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Figure III-99 

Turn Phase Run C Lateral Distance DY 
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Turn Phase Run C Yaw Difference 
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Figure III-101 

Turn Phase Run C Longitudinal Position DX 



5 ♦ Save E ffects on Velocity C ontrol 



The final testing- procedure involves validation of • the 
speed control system in the presence of waves- This 
perturbation testing continues that started in section A. of 
this chapter for heading control. In chapter IX the wr 

force was modeled thru the intermediate force IF32' asr 

% 

IF 2 2 = KC1 •D2+NC2+KC1*fiX 

By introducing the force in this way, a severe limitation is 
placed on the magnitude of the force. In the mariner model 
used,, the KC1 coefficient (XDELfi) is considered negligable 
or,, at best, only 0.00005. This translates, in the original 
equations of motion, to a maximum speed perturbation of only 
0.C355 Jets, for the wave amplitude chosen. The second 
drawback of this method, with even greater consequences, is 
that the perturbation is introduced before the control loop. 
Delay of the wave perturbation is produced making it out of 
phase with the other wave force {WY) and moment (WN) . 



WX 




Figure III-102 

Block Diagram of Wave Introduction in Speed Loop 



In order to bring about uniform introduction of this 
wave force, its effect is inserted just past the integrator 
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of the speed control loop as shown in figure III-102. This 
is coded in the DSL simulation program as: 

CD0T2 = INTGRL (U02,SPDEBE*LUC)+KS2*WX 



A value of -1.0 for KS2 will give a maximum wave 
perturbation of 0.85275 kts. (a much more realistic 
perturbation for the high sea state simulated) . figure 
III-103 portrays the speed desired and speed acquired for 
the approach phase in the presence of sea state. From this 
it can be seen that the speed acquired is very dependent 
upon the sea state present. The control law, however, 
presents a very stable reference for the speed loop which 
gives an approach longitudinal position (DX) plot 
indistinguishable from that of figure III-81. More 
prominent perturbation results are evident in the turn phase 
plots of figures III- 104 and III-105. The speed response of 
figure III-104 allows a maximum longitudinal position 
excursion of 9.5 feet (0.018 normalized position), as 
compared with 8.286 feet (0.0157 normalized) in calm sea. 

These results show that the speed control system is 
very stable and corrects well for large external 
perturbations. 



195 




196 



Figure III-103 

Approach Phase Speed Response in Waves 
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Figure III-105 

Turn Phase Longitudinal Position DX in Waves 



IV. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLCSICNS 

The results of this design study have been most 
gratifying. The basic concepts initially perceived for the 
HAS control have been realized. The decoupled ship control 
in the BAS environment is a viable and plausible idea. This 
thesis contains a workable system for implementation of 
computer controlled RAS. The achievement of 2.3 foot 
maximum excursion for lateral distance while both ships are 
in a turn, and longitudinally offset by 53 feet is a 
phenomenal achievement. Having this kind of accuracy in BAS 
operations, can vastly increase the safety of this 
complicated and dangerous maneuver. 

The approach phase of HAS can be a very hair raising 
experience. Night replenishment and sea state complicate 
the ’’seaman's eye" method now employed in the fleet. Having 
a system that automatically handles the approach regardless 
of the adverse conditions can, again, do nothing more than 
increase the safety of the RAS maneuver. 

Schemes for computer control of nonlinear systems and 
the purposeful introduction of nonlinear control laws are 
becoming mere practical with the technological advances in 
micro processors. The ever increasing number of U.S. Navy 
ships with computer systems installed, makes digital 
computer ship control realizable in the present time frame. 
A good micro computer or an existing installed computer 
{such as one used for the NTDS system) can be used in this 
vein. Procurement of the hardware required for this RAS 
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system can be dissipated over time periods contingent cn the 
funding available. The supply ship requires only two 
reflectors fcr the range and bearing devices stationed on 
the receiving ship. All ships can be outfitted with such 
reflectors at a minimal cost, while the bulk of the hardware 
can be introduced to the ships at regularly scheduled yard 
periods. 

In the initial conception of this thesis, a section was 
planned fcr open ocean maneuvering, after seme research on 
this facet' of ship control, it was determined that work in 
this area has already been documented[ 24 ][ 2 5 ]. The 
existence of NTDS outputs for station attainment and single 
ship control systems, made design in this area redundant. 

The concept of integrated centralized ship control has 
been in the background for over a decade[ 26 ]. Although 
given a lew priority due to funding considerations, its 
implementation seems to be just around the corner[ 27 ]. 
However, a review of ref. 27 indicates a lack of HAS 
capability. Whether this is an oversight in the article or 
neglected in the design criteria is unknown. If it has been 
neglected in the design, a very real problem has been 
overlooked. The recent incidences of ship collisions while 
conducting HAS [ 28 ] emphasizes the need for inclusion of 
this very dangerous maneuver in the "Integrated Eridge 
System." Lack of technology can no longer be used as an 
excuse. This thesis and other research reports [ 29 ] have 
advanced the implementation feasibility to a level that 
cannot be ignored. With these projects finalized in 
practical terms, their incorporation into fleet use is the 
next imperative step. 

A major effort in this area must be made. The ever 
increasing complexities of today's naval ships and the leads 
being placed cn the officers and men are such that computer 
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control must be used; and used now! We cannot afford the 
luxury of time to prove these sy stemi^Xcrth , but must make 
concerted efforts to get them implemented before the lives 
of 300+ men are lost. 

Whenever a complicated system such as a ship in the RAS 
situation is encountered, many facets have to be 
concurrently analyzed. This fact has caused inclusion of 
many diagrams in this thesis to illustrate the total 
picture. Each run, with a different condition, reguires 
many plots to analyze the differences in the responses and 
the causes of the differences. The computer programs shown 
dc not reflect the actual run times in the JCL shewn. As 
many as twenty plots were output in these programs in the 
times listed. Analysis of the actual computation times show 
that the algorithms run considerably under the time reguired 
for real time operation. The sampled data rate used in the 
simulations was 0.11 seconds. This is well within the 
realizable data rates available in even the slowest of 
today's computers and microprocessors. The thrust of this 
consideration is that there are nc problems envisioned in 
converting BAS simulation to real world RAS control. 

E. BECCK HE HE ATICNS 

In the beading control design section of this thesis, 
the need for a completely adaptive gain scheme was cited. 
Again in the velocity control section, when a longitudinal 
offset was introduced, this need became even more evident. 
The first and most important recommendation for further 
study is the development of just such an adaptive gain 
scheme. 

The linear equations of motion should be replaced with 
nonlinear eguations to validate the control designs advanced 
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Along with this. 



in this thesis. Along with this, the hydrodynamic 
coefficients for the Navy's modern ships are required to be 
able to design these control systems for todays vessels. 



It is further recommended that a concerted effort be 
made to obtain data on the interactive forces and moments 
between ships of dissimilar types and sizes. These forces 
and moments must also be available for sea state conditions. 
In fact, the whole area on sea state effects cn the various 
ship types in the HAS situation and in open ocean 
maneuvering needs attention. Not enough data was available 
for this researcher to be able to pinpoint sea state effects 
on ship hulls. Since replenishment at sea is rarely 
conducted in the sterile condition of calm sea, these 
considerations are of utmost importance to allow testing of 
any control system in the s-imulation stage of development. 



The intent here is not to imply that the control 
systems portrayed in this thesis are the best for the HAS 
scenario, but that the procedures used can be applied to any 
control scheme desired and bench marked to the ones 
contained here. As previously mentioned, much meaningful 
research and design must be accomplished to allow system 
reliability and, mere important, system acceptability by the 
officers and men who will ultimately trust their lives to 
it. This is a task that must not be taken lightly. 
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APPENDIX A 



Due to the lengthy nature of the computer programs 
presented in this thesis, many functions and subroutines 
were developed tc simplify their presentation. This 
appendix lists these functions and * subroutines in 
alphabetical order. The computer programs reference this 
appendix and indicate the placement of the reguired 
functions and subroutines. 



A brief description for each listi 
the reader in determining their purpos 
following is a listing of the functio 
contained in this appendix in the order pr 



ng is given to aid 
e and use. The 
ns and subroutines 
esented: 



S CE ECCTI NE BOXPLX 
FUNCTION DEGRAD 
FUNCTION DELAY 
FUNCTION FE - RUN A (FEA) 

FUNCTION FE - RON E (FEB) 

FUNCTION FE - RUN C (FEC) 

SCEECUTINE HDGRAS 
FUNCTION KE 

MAIN PROGRAM FOR FUNCTION MINIMIZATIONS {MINIEXPX) 



SUEHOUTINE RBMEAS 


FUNCTION 


RKLDEQ 


SUEROUTINE SLOPES 


FUNCTION 


SPINIT 


FUNCTION 


SPDCTR 


FUNCTION 


SPDOFC 


FUNCTION 


SPDREC 


FUNCTION 


SWCL 


SUEROUTINE SWTCH 
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SO EBGUTINE SWTCHF 
SOEECUTINE TRANS 
FUNCTION XLiaiT 



2 04 



SUBROUTINE EOXPLX 



This subroutine was used for all optimization runs in 
heading ccntrcl and speed control. It was programmed 
locally ard is part of the IBM 360 SSP library at the Naval 
Postgraduate Schocl. A full explaination and description is 
shewn in the first few pages of the subroutine listing. 
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FUNCTION DEGRAD 



This function was programmed to convert degrees to 
radians and radians to degrees. A third purpose is to 
convert degrees to a range of C-360. It is used extensively 
throughout tie programs, functions, and subroutines listed 
in this thesis. 
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FUNCTION DEGRAD ( NDRFLG, NSHFLG t FUNCT ) 
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FUNCTION DELAY 



This function is used as the time delay in the speed 
ccntrcl optimization runs. It was designed to be used as an 
equivalence to DELY in DSL simulation. The following 

variables are defined: 

E(I) is the storage array (should be initialized 
before the first function call) 

K is the delay step count 

SPDDES is the variable to be delayed for K steps 

P is the flag for delay or no delay 

P>0.0 delay SPDDES 

P<0.0 function output equal to SPDDES 

The function stores the input value (SPDDES) in E ( M) and 

decrements the value in array E (I) at each call of the 
function until the value is in the position of E ( 1 ) . The 
value is then output from the function delayed K intervals. 
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FUNCTION DELAY 
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FD NOTION FE 
FEA 



RON A 



This function is the simulation for heading control 
optimization of the approach phase. It is called by 
subroutine EOXELX. The integration step size is 0.04 with a 
final tine cf 20.0. In this function all initial conditions 
are set to zero except initial geographic location and 
speed. The reference ship maintains a straight course and 
the control ship starts its approach 5 ship lengths astern 
and 0.4 ship lengths laterally displaced to starboard cf the 
reference ship. 

The function is referred to as function FEA in the 

text. 
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FUNCTION FE - RUN 
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FUNCTION FE - HUN B 
FEB 



This function is the simulation for heading control 
optimization of the turn phase. It is called by subroutine 
BOXPLX. Ihe integration step size is 0.04 with a final time 
of 20.0. In this function, the following initial conditions 
are ncn-zerc: 

control ship rudder angle D2D S I (14) = 8.7 degrees 

lateral displacement Y ( 1 1 ) = 0.2 
reference ship's speed U01 & Y{4) =1.0 
control ship's speed U02 & CDCT2 = 1.5 (after first 
step becomes 1.0) 

The reference ship's rudder is activated tc 5.0 degrees 
between time 4.0 and 5.0. The runs were for port side 
replenishment . 

The function is referred to as function FEE in the 

text . 



230 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

ro tn so oo O' o m 4 - in nO r*- oo o ^ ro in >o co on o m oo O' o m 

^^^^^ l -i^^*^ f «,cNj(\jcNjc\in4c\ifsjc\if\irvJrnrornrncnrn(r»rnrnfn4*4 ,, 4“ 4* 

COCDCOCDCOCDCOCOCDCO CDCDCDCOCDCDCQ CD CO CD CO CO CO CD CO CD (DCOCOCOCDCDCQCOCOCOCOCOCOCDCOCOCO 
rsl rsj M *4 M r*4 rxj fsj M rvl rxj rM fsj f4 pv4 rxl M i >4 isj rxj rxj rsl rxj M rsl rxj M i>4 rvl rxi ivj rxj fvl rxj rx4 rxi rsj rxj rsj rM rsj i>4 rxj 
UJ LU LU LU Ui UJ UJ UJ UJ LU LU LU LU UJ LU LU UJ UJ LU LU LLi LLJ LU ULJ UJ LU LU LU UJ LLI LL1 UJ LU LU LLJ LU LU LU LU LLI UJ LU UJ 
lLlLU-U.lLLLLLUaLLLU.U-U.LLU_U.U-U.U.lLU.LLlLtLULULlLLLLLIJLtUU.LLU.U.LLLLLLlLLLtLULLL 



O 

QC 

•* 

CO 

z 

UJ 

CO 



CO 

z 

UJ~* 

COO 

a: cm 
CLt- 

oo 

a 

2X 

a 



oo 

2 CM 





d— 










O 












ULX 










z 






















0 












< - 










0 












«•*» 










LL 












COO 






















<L CM 










LU 












w 










CD 












Z ►— 
















Z 






OO 










O 






O 






*-*Q CO 










1 — 






*— 4 






»-> h- 
















K- 






o z 










CO 






< 




o 


Z •* UJ 










Z 






O 






o— w 










»— • 






O 


CO 


z 


LLOOO 










< 






O 


z 


o 


CMf-*— 




CM 






O 








0 


Qd 


h— - -LL 












CO 




O 


1 — • 




C 0 >QLL 




< 




2 : 


LL 


Z 




w 


L- 


1 


O 1 — UJ 








LU 


O 


O 




X 


< 




— O -QO 






in 


V- 




*— ♦ 




CL 


rxj 


UJ 


r *4 — • -O 


r- 43 


<M 


43 


O 


Z 


►— 




< 


M 


LL 


w LL 00 X 


mN a>- 4 oor- cm 


< 


I s — 


O 


O — 


t— » 




a: 


O 




Ow ^o 


4 - r\i in sO cm tn cm cm I 


4 - 


M 


t- 4 -%CM 


Q 




0 


< 


z 


UJ psJh- winrsJOinOrOOcMCO^cMO^cM 


COOOZ 


h— — *w 


-Z 


4 


0 


OCMI^-h-i 


o 


U-Z OS-*hOOOOOOOOOOO<M 


• • 


•O 


<twivi— • 4 "Q 


— »o 


LU 0 0 


• • *H- 



QZO <00000000000 ♦ 0 < 00 (NJNKCN 0 N 
h- ZmQXZ •*•••••••• »0 •CM O • • » ll '-'ixa 

o a v-hoo >00000000000 i o~«o imioocou.ii >4 _j 

z — <iuo* a — * ii *-* co ii < 

O KDZ O II li I* II II II II II II II II II II < II St II II II h-<S)Z ll N* 

UL OOLU-J0C UQLLI ZZUJ Of- 

•-*CMrMPvJCMfncO<CDO II UDKKOQQUiUJoozaQM 

D>*-iuj>^^rMcM^r-«cM<Mrricni£*:^ '-‘cmQuoh-oll^ 

lLUjQctX<CD<CC<CD<CD<CDXXXQXXQQXQQMa:2a:>»H 



•— »o 

II 

o II 



Ow 

o> 



O • *OO00m 

oo z 

—I II ll II ** 
< II II 

h--*— -O— «4UJ 
•— 1 m \Q ^ JZ 

~>»>>a 



o o 



o oo oo 



o 



231 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

a^o^osjrn^m vor^oo o'O— ♦cMcovt-in-or'-co c' 0 ^-«(\jmvfm«or— ooa'O—* 

4 - ^ ^ ^ in tn in in in in in in in so nO >0 vO vO ^0 vO vO >o sO ^ ^ r^- oo co oo oo co co go oo qo oo cn 

CO CO CO CO CO CO CO CD CO CO CO CO CO CO CD CO CO CO CO CO CO CO CO CO CO CO CO CO CO CO CO CO CO CO CO CD CO CO CO CD CO CO CO CO CO CD CO 
h 4 ivjfM rvj rvi rvj rvi rsl hj r ^4 M pM fvl r^i rxj M r^J iv* rvi M M ivi fvi ivj r^j ivj rvi rvi rvi hsi r^j m M r*u rvi rsj M ivi r^j tvj fvi r^j ivl Psl *nj rsJ rsj 
UJ LU UJ UJ UJ LL1 UJ UJ UJ UJ UJ liJ UJ UJ UJ LLJ UJ UJ UJ UJ UJ LLI UJ UJ UJ LLJ LLJ UJ LLJ LLJ UJ UJ UJ UJ UJ LLJ UJ UJ LLJ LU LU UJ UJ UJ UJ LLJ LLJ UJ 
l^U»U-LLtX»lLLLU»lLU»LLLLll.U. LLLLLLLLLLILILLL LLU_U.U-LLLLli-LLLLU-LLLLLLLLLLU.lLLLLJLLL LLLLLLLLU.IL 













• 
























in 










— 










— 




11 










r\iosi 










cm 




00 










<W'«^ 














UJ 


— 








>> 










> 




O 


1 — ^ 


















#• 






UJ 








ZOO 










>--4 




Q — 


H- 








H-.O 










oz 




00 


-J 




Q 


0 


000 










<>• 




— UJ 


O 




~HCM X 














*• ^ 




—0 


•» 




LLLL — 




— — 










xcm 




0 -* 


X 




*-44—4 —4 


CM 


— •—4 










Q>- 




• Q 


UJ 




XX c\j 














O 




in •* 


0 




+ + <r 


<1 


>• >• 










•» •» 




•0 






# 




1 > 










o*—* 




UJ • 


QJ 




mm ^ 


—4 


— 










>> 




-JO 


•*X 




CM 


p4 


— ^ 








o 


Q>- 




• m 


SI* 






M 


(MCMcvj 






uo 


t— « 


•* •» 




H- 


UJ»-H 




•w* mx 


X 






«—» 


2 




o> 




—O 


H-UJ 


— 


-4CMLL 1 


1 


>» 




Q 


a 


< 


xo 




• • 


-JCQ 


Q 


C\JCMi— iCVJ 


>— 4 


•» 




(\J 


M-l 


fxj 


a< 




OO 


OH 


1 -^ 


UCOCOXCM 


^♦rnoozo 


C\J 


a 


H- 


»—♦ 


#• •» 




ZroQ 


1 -J 


O 


Zl l+< 


<mO M •• 


> 


•» 


<t 


— 1 


— x 




<C 1 — i 


i«~0 


•* 


x-**— «. 


& <L UoOO 






rvl 


<inocMQuj 




• wQ|_^ 




+ 1 <■■■< 








N-4 


N4' 


— cr 




— t— 1 


— »C£) 




•—4 


,OJ 


ro 


C\J 


—4 










o*-**-< 


5ZX 




0 




<1 


— 


1 1 ^ 




• SOOnhX 


-ft S* # ■* # X*- 


»x 


— w^qc: 


* 


Q 


M 


inZ 


— • 00 






— Q<IcOO'^njrO^> 




x>>o 


<t 


<t 


h-O 


• »— « 


O-*00UJUJ 




• •-J-JXII CM<XXX— «*^m 




UJ 


X 


lOOOON • 


om— ♦ 


»^Q.lvJ 




OUJXQ 


— *QCXXXcQCQCQ II II 


II 


II II II 0 


X 


DO • 


♦ . * 2 :^. 


• • UJ 


-w<Om 


0 


O 


II — 


wO 


1 1 1 









• • »UJOOOOm II O»-0C—l-JO • II • II II <M>-UJ || || || «l II 

HHOOQ II z >— oO< •Oa< H- — • Q II II II — ♦(NJrOvfinsO 02 

r-* II II || II Q II II <M< II II HO II —♦—«—« CM 

ii ii ii i) ii uj— h-h- jjk it n a.ujwooujcoH u ^cxim^— ♦— » h-i— i— i— h-h-Qcl^ 

II — IC\J— ♦C\JLUv|"^(\jO^ 0 OO— I— II ►— »Q hhhOQ U.lLU-HN^OOOOOa 3 Mii. 

OooO-<rg»2ZQ.'-OOO w >X«2 hhX^U iJjQhhmmmmmhOOQOQO<Ih 

Za:wQ 0 Q>»>^>D 3 oQ 0 QUU M l- 0 “ 3 ^Q M t 



<M 



u u o o 



u 



232 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

cnj m tn >o oo o — 4 m in >o oo o ^rsj m in vO oo o o rvj i'o vj* in vO oo O' o ^ c\j ro ^ m >o r- oo o 

^O'O'^O'O'O'O'OOOOOOOOOO^- i^^^^^^^CNjcvjc\J<\JcucNjcgc\jc\jf\irornrnrornrornrorOcn 



CDCOCDCDCOCOCOCOCOCOCOCOCOCOCOCQCOCOCOCQCOCOCOCOCOCOCOCOCOCO ^rnr prf^ rnrr^rOr p rOrfirn^O ^ ^H > ) r r i rnrn 

rsJisjisjrsjfsjfNjiNJ^^^^hsJiNjfNj^rsifv4^rsJrsiisjrsjisjrsJrvirNjf^jrsJiN4rsj^ 

LLIULI 1 U* U * | j i |j 1 1 iimiJilJiiiLUiiiLULUilJUJilJLULULUuilUii in iLUmLUMjnnim pi Minmi j n iujii | i|ju ju_i ijji i MIJLU 11 Mil 

U.U.U.U.U.li.U.li.LLU.ULULULLLU.LLLLLLU.LLU.LLLLU.LLLLLLLLLLUuli.U.LLLLU.li.LLLLULlJLlLU.LLLLli.LLLLLL 

•ft 4—4 

CM 00 

CD Q. 



(V 

cc 



O 

O 

< 



CO 00 
CD Q- 



CO LL 



O 

LL 

O 



O 

cl 



oo 

Q- 



00 

> 



(MCM Q 

-4 CM ^ 
ILIi- — 

►-44-* 4^ 

XX — CM 

+ + <M < 

— o * 

On CTO <M 
—CM CM 
>>-*- 4fO«-4 
# * OULX 

— iOJZJh-. | 
(MCNJ -XfM 
OJOCOcOX + CM 
ZZ I I Qrg< 
>X— <1— -* 
+ 4 -NNwQN 



00 CD 




OsJ 

Z 

> 



00 



00 

z 

UJ 

oo 

a: 

«• 

r\J 

CD 

CD — 

•O 



— Z > 
» •» 
O - — 

<<M sO 
->- — 
x> >- 

O •» •* 

•.>— *in 
^>-co— 
Q •*— ’> 
-» - 
XQ — 
Q<OCM 



CD 

OS 

•r? 

OO — 
nOQZ 



-4>>- 

X—— 

» OOZ 
»-<O m 

*-4000 

<£#> # — •* 

— XO>- 

ojQ — - 

cm cm — ql r> — cm O O in «d — <c q ^ 
OQ»hO M O w QQ^>-> , -<^ 
* # # #Ott X— XUO>>— QC 
cOO^NQiO^yw | | 00 O 00 

XX-«-4a.m ii it — voujllki 

XXCQcQOOCQ II II II O^ZOuQUJ 

| » | — »^- 4<0 S 

II If II — — .o-*— cx: — I II CD 

II II iH*-ooon-4~4>->i— oo cl 

CM CM CM CM 

CMrOCMCMh- CMh- I— •— I — ♦— II It. 



CMQC: 


m + o 


— * 












CD - 


l oa: 


UJ 






A 






-o 


co • 


t— 






h- 






rvjQ 


ZOCO 


-J 












a: - 


ooz 


o 






— 






-v> 


omo 


v» 






O 






-»Z 


COO II o 


X 






1 — 






COUJ 


LL II CO a 


UJ 






o 






c000-* 


1 — CDZ — 


1-0 












— 


OSDO 


-O 






o 






-tZSO' 


ooo • 


Q-J 




— ► 


►— 






CO — ^ 


COQQO 


••x 




— 








~o> 


+ 0 rO 


X* 




-5 


1 — 






—•CM -OO — * - 


LLIUJ 


— 


<— 


o 


— 




QdcDi- 4 (MlU <-»00 


H-cO 


o— 


h- 


o 


—• 




— — » 1 1 ■ Q O • • 


•Jh- 


<M>- 


O 


X 






(/)OOQw *00 


0—1 


QQ 


— o 


•» 


— 




4 — • —4 — O 00 O CO CO 


i a 


-<tUJ 




X 


—X 




-co o co cl oo —• i a 


— & 


—* 1 1 — 




•ft 


- 5 — 




2 »<•» 


| —♦ | ^f\Jh -0 


-Ooo 


— UJ 


— xi- 


— — J 




'-'OQCOCM • *t— O 


—4^: 


hQm 




H— *— 4 


xo 




^OcOOk t~— 1 


XX 


— wO 


— CD 


— O — 


— z 




00 <UJLLXO-JZ 00 


IHH 


-•Ooo — * 


^•UJO 


i— or 


vt — |L0 


o 


C -ox • • 4 -»h--J II Nt<cO II ro 


*- 4 — J 


— UJUJ 


-40 


t- 


<*< 


>C 0 C 0 -I-IX - 


-»a:< 


—CO II 


UJ— JO 


-Z 11 





OO II II SSXQ — 

Cj - ll or? II >*->UJ II rO> 

Xh-Qco oo ii ii — • a — i 

3 SCMU.COQO UJ— II UJ— II 
rOOUJCQI— III 



U.LL^(MO0OOOOOO jJ2J(MJQOOr hf-HOO oOO^ 

MM MMQMQQQOQX>««CO<UJOQOlLliJJ-iOlM(MMQcD 

xxxxox>»»oquu>ucqoococqQhmoqo>ooo>o 



ft-lQ 


— ICO— 1 


-4 tO 


— 1 


II — 


OQX 


II — 


o 


-5 II "5 


O CL 


— 5 II — > 


z 


cO— 


II 

It II 




00 


OO 

0—0 


o 

Oh- oo 


-50 

0—0 


II 


oxx 


1— Qfsj 


o>> 


h— 



m 



o o o 



o 



233 



ooooo oooooo oo o 
O f^*oj cn >4* in co & o *-« nJ co 

>t ^ ^ >s- ^ >4- in in in in 

OOOOCOCOOCOcOCOcOCOCOCO 
rsi rsl in* rvl tvl rvl inJ t-M rsi r>i rsj fM rsj ivl 
wmmmmmmmmujmmmm 
U. UL U. U. Or UL IL U. U. U. u. U. UL U. 



I! 

—X 

- a 

s< 

m 

o- — * 

CO **C0 

a qo • 

c£ # m 

O-in^ 
— «UL 
ZUL - 

Q — 

H-»* II 

hlio 

<><M 

c^QO 

o<x * 

LU 

f— • * 

z ~x> 

M GO * 

•in 





Q 


•►in— « 




CM 


O^U- 




a 


•LL * 




•» 


— ♦ •• 




CM 


• II 




O 


1 — 11 CM 




3 


030 




< 


C03 




> 


mO< 




•» 


3 ->» 






-J-— - 




a 






2 


>— » * 




<t 


LLiCO 




>► 


OWL • 




•» 


m in 




X 


ZZ-* 




a 


a:Om 




< 


3— - 


o 


•» 


h-H— 


h- 


> 


mo il 




Q 


q^Zh 


a 


<t 


30 


sO o 


•• 


Cj*UL2 


•» 


-? 


LU <1 


-*<\l -* 


CO 


Oh> 


o * o 


a 




h* in • 




XX •* 


a — o 


— 


c^mx 


'——.CO CM 


O' 


in 


O • *. • 


•» 


W V •» 


0»— •'O 


-5v0 




Z | OOJCG— 


•• 



LOoO • a z i— h>co 

rsim i—o met << • 
it O.^H II i— 3 zzin 
hh o i— « i— a:a^a 
hU. 0 CHU.OmcXUJ OOa .2 
OH-<3to*-»OLL3ct mm m 



in >o fn* 



GO O' 



o 



234 



FUNCTION FE 
FEC 



RUN C 



This function is the simplified simulation for speed 
control optimization of the switching function SW. It is 
called ty subroutine BOXPLX. The function shown is for 
approach speed of 1.1 and a replenishment speed of 1.0. The 
runs were made for various realistic combinations to obtain 
an optimum switching curve. 

The rut used a sfep size of 0.04 and a final time of 
10.0. The two ships were run linearily with only the 
longitudinal direction and motion of any concern. 

This function is referred to as function FEC in the 

test. 
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SUBROUTINE HDGRAS 



This suhrcutine was programmed to calculate the desired 
heading (ESIDES) for RAS heading control. It uses the 
outputs of subroutine RBHEAS to calculate this heading with 
gains RSENS and WTSENS. The large number of outputs in the 
subroutine call statement were made for ease of DSL printed 
output fcr tracking of simulation accuracy. 

The subroutine also incorporates a loop to avoid 
computer precision problems in the ARSIN function. 
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SUBROUTINE HDGRAS 
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FUNCTION KE 



This function is required by all optimization runs. It 
is the function that contains constraints for subroutine 
EOXPLX. No constraints are present, consequently function 
KE is set egual to 0. 
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FUNCTION KE 
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BAIN PEOGBAM FOE FDNCTIO N MINIBIZATIONS 

MINIBXPX 



Ihis is a generalized program which calls subroutine 
EOXPLX. Its main purpose is input and output of the values 
reguired in the optimization runs. This is referred to as 
MINIEXPX in the text. 
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MAIN PROGRAM FOR FUNCTION MINIMIZATIONS 
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SUBROUTINE RBMEAS 



This subroutine measures the range and bearing of the 
forward and after stations which is required of subroutine 
HDGBAS. This is done with trigncmetric functions as shown 
in chapter II. The subroutine is specifically designed to 
circumvent any ambiguities usually associated with these 
functions. 

It is the basis of the decoupling of the two RAS ships 
that this thesis is based. 
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FUNCTION RK1DEQ 



This function is the Bunge-Kutta-Gill forth-order 
integration used in ail optimization runs. It is programmed 
locally and is part of the IBM 360 SSP library. A full 
explaination and description is shown in the first few pages 
of the function listing. 
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SUBROUTINE SLOPES 



This subroutine contains the table look-up and 
interpolation scheme for the interactive forces and moments 
presented in the HAS environment. It is long and must be 
pre-compiled for most of the DSL simulation programs shown 
in this thesis.’ 
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FUNCTION SPXNIT 



This function was designed to aid in initialization 
problems associated with the DSL function DELI. The effect 
is that the function initializes the delay loop until it can 
be self-sup pcrtive. 
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FUNCTION SPINIT 
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FUNCTION SPDCTR 



This function is the speed control used during heading 
control development. It is used directly as the speed of 
the control ship with information of the speed of the two 
ships and the longitudinal position ADX. It contains a 
linear function at ±1.0 ship lengths to a dead zone of 
±0.001 centered about the alongside position (0.0) . 
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FUNCTION SPDCTR 
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FUNCTION SPDOFC 



This function is identical to SPDREC except that 
ability tc offset the alongside position (O.C) 
incorporated. This is the speed control function in 
final develcpaent form. 
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FUNCTION SPDOFC 
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FUNCTION SPDREC 



This function is similar to SPDCTR except that a 
switching function is incorporated. This is the function 
used for optimization of the switching function and is used 
in the velocity loop simulated in the velocity control 
section cf chapter III. 
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FUNCTION SPDREC 
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FUNCTION SW CL 



This function contains the fifth order polynomial curve 
fit fcr the optimal switching position of the speed control 
lcop. Its range of values for SPDDIF are 0.1 to 1.0 
normalized speed difference between the two ships. 
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SUBROUTINE SWTCH 



This subroutine contains the gains and mechanisms 
reguired for the adaptive gain schedule developed in this 
thesis. It includes the optimal gains obtained from the 
heading control optimization runs. 
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SUBROUTINE SWTCH 
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SUBHOOTINE SWTCHF 



This subroutine is identical to SWTCH except that the 
turn phase gain 7PBG is relaxed to allow for offset 
longitudinal position placement. This is the adaptive gain 
schedule in its final form. 
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SUBROUTINE SWTCHF 
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SUBROUTINE TRANS 



This scbroutine takes the lateral and longitudinal 
geographic displacements and converts them to actual 
displacements referenced to the control ship r s head. This 
is done tc gain a more realistic reference for subroutine 
BENIAS, subroutine SLOPES, function SPDCTR, function SEDCFC, 
and function SPDREC. 
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SUBROUTINE TRANS 
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FUNCTION XII MIT 



This function was developed to allow the LIMIT function 
of DSL tc be incorporated in the optimization runs. It is a 
saturation amplifier with a gain of 1.0, and upper limit of 
UL, and a lower limit of XLL. 
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APPENDIX B 



The final fcrm of the simulation program, with all its 
subroutines and functions, is a very complex and complicated 
maze. Tc aid in following its progression, this appendix 
contains a detailed block diagram of the program with each 
variable listed in its computer variable name. Each page 
contains a functional part of the simulation with inputs and 
outputs shewn cross referenced to their origin and 
destination . 

The following is a list of the block diagrams contained 
in this appendix in the order in which they appear: 

Ship A (Reference Ship) Simulation 
Ship A Heading Simulation 
Ship A Speed Simulation 
Ship B (Control Ship) Simulation 
Ship B Heading Simulation 
Ship B Speed Simulation 
Subroutine RBMEAS 

Bange Measurement 
Eearing Measurement 
Subroutine BDGBAS 

Heading Control Loop 
Auxiliary Functions 
law Conversion 
Coordinate Conversion 
Feedback Loop 
Rudder Modeling 

Ship A Rudder 
Ship B Rudder 
Have Generator 
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Wave Direction 
BaEdcm Variable 
Wave Encounter 
Wave Components 
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Ship A (Reference Ship) Simulation 








fjr SWAY2 
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Ship B (Control Ship) Simulation 
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Subroutine RBMEAS 
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Subroutine HDGRAS 
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Auxiliary Functions 
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Rudder Modeling 
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Wave Encounter 






APPENDIX C 



Throughout this thesis subroutine SLOPES has been used 
to output the interactive forces and moments between ships 
in the RAS situation. This subroutine, adapted from ref. 
11, does net contain a complete picture of the circumstances 
envisioned. In particular, ship’s speeds ether than the 15 
kt. operating point and different ship lengths are not 
accounted for. 

As stated in chapter II, the speed modification factor 
can easilj he applied for both ships at the same speed and 
other than 15 kts. with the following expression: 

z 

SPDP = CDOT 

Ships replenishing with different lengths can also be 
incorporated as shewn in ref. 1. 

Subroutine FAMIC listed in this appendix incorporates 
these twe ideas along with a better method of deteraining 
the interactive forces and moments. The curves of figures 
11-11 and 11-12 were quantized every 50 feet of DX for all 
the DY curves shown. These points were then used in the 
NEGS XDS-S30C digital computer and AGT-10 graphics terminal 
tc obtain a family of best fit curves. The test fit 
criteria is based on the sum of the error squared at each 
quantized point (modified somewhat by this researcher's 
evaluation of test fit between points to eliminate spikes 
and ether anexalies) . The results of this curve fit process 
is summarized in tables C-1 and C-2, which includes 
tabulation cf the best fit criteria. These polynomial 
coefficients are based on the DX distance and are coded in 
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Interactive Curve Fit Polynomial Coefficients 
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0.176 




3.96169 


0.252 



Table C -2 

Interactive Curve Fit Error Analysis 



303 



subroutine F A KIC as YY(1) thru YY (11) and YN(1) thru YH(11). 
An interpolation algorithm is used to determine the forces 
and moments at DY points between the curves of each family. 
Although all the computations are based on measurements from 
the control ship (ship #2) , the interactive forces and 
moments are also computed for the reference ship (ship #1) . 

Figures C-1 and C-2 are the interactive forces and 
moments output to show comparison to figures 11-11 and 
11-12. The speed of this run was the operating point of 15 
kts. The ships are of equal length (527.8 feet) . 

Linear interpolation of the interactive curves for 
greater than 150 feet DY distance is accomplished from this 
150 foot curve to a value of 0.0 at 200 feet. It therefore 
assumes no force and moment are present outside the 200 foot 
range. All forces and moments for DY distance of less than 
50 feet are taken as that of the 50 foct curve. These two 
endpoints are by no means exact, but will suffice until more 
detailed data can be gathered. Another inexact endpoint is 
produced at the curve families limits of ±550 feet. At 
these points, the forces and moments are forced to 0.0 since 
detailed data outside of these limits was not available. A 
side effect cf this abrupt truncation will manifest itself 
in the instantaneous commencement of the forces and moments 
during the approach phase run. The endpoint variations in 
seme of the curves of figures C-1 and C-2 are due in part to 
the curve fitting routine used, but mostly to the 
differences in computer precision. (curve fits were 
calculated on a 11 digit precision XDS 9300 while the curves 
were plotted on single precision 7 digit IBM 360/67) 

As previously mentioned, the speed modification for 
other than the operating point cf 15 kts. is only completely 
valid for the situation where both ships are at the same 
speeds. Since this thesis considered an approach phase 
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Figure C-l 

Curve Fitted Interactive Y Forces 
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Figure C-2 

Curve Fitted Interactive N Moments 
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where the control ship enters the interactive field at a 
speed guite different than the reference ship, some 
modification of the interactive effects should be 
considered. However, exact relationships are not available 
to compute the required modification factors. 

To dispel any problems with the design of the heading 
control system, the worst case speed modification factor was 
chosen. This factor, in effect, considers that the 
interactive forces and moments are derived from the control 
ship. This is accomplished in subroutine FAMIC with the 
following fcrtran expression: 

SPEE 2 = CD0T2**2 



As stated in chapter II, it is felt that it is more 
accurate to consider the interactive forces and moments to 
be modified by the speed of the reference ship, and can be 
coded in subroutine FAHIC as: 

SE IE 2 = CDOT 1**2 

fiith the scenario followed throughout this thesis, this 
expression would equate to unity throughout the HAS 
situation, since the reference ship is maintained at 1.0 
normalized speed (15 kts.). 

For the sake cf error analysis, simulation of the worst 
case modification is performed. This gives rise to forces 
and moments 2.25 times what they were in the rest of this 
thesis during a portion of the approach phase when the 
normalized speed of the control ship is 1.5. Figures 03 
and 04 show the interactive forces and moments for the 
approach phase cf the simulation. The comparison plots 
which appear in chapter III as figures III-24 and III-25 
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Approach Phase Curve Fitted Y Forces 
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Figure C-4 

Approach Phase Curve Fitted N Moments 
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Approach Phase Geographical Plot From Modified Interactive Effects 




illustrate tbe extent of the changes. Most notable is the 
smoother output of subroutine FAMIC. This mere 
realistically portrays the interactive effects in the HAS 
environment. Figure C-5 portrays the geographical plot 
which compares with figure III-26 without speed 
modification. Although differences exist, figure C-5 
illustrates that the interactive effects speed modification 
factor for the worst case does not drastically alter the 
approach phase outcome. The heading control system design 
is still valid in the face of these changes. 

For reference, figures C-6 and C-7 show the interactive 
forces and moments in the turn phase as calculated by 
subroutine IAMIC. Figure C-8 is the turn phase lateral 
distance plot produced. It can be seen from this 
illustration that the maximum excursion is 0.0056 normalized 
distance (2.56 feet), well within acceptable limits. 

In summary, the designed control system will accomodate 
even the worst case modification of the interactive effects. 
This insensitivity to a large range of perturbations, makes 
this control system a more viable design for actual ship 
ins tallaticn . 
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Turn Phase Curve Fitted Y Forces 
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Figure C-7 

Turn Phase Curve Fitted N Moments 
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SUBROUTINE FAMIC 
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POLYNOMIAL POWER PRE-COMPUT A T I ON FAMI 
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COMPUTES PROGRAM #1 



This program incorporates the ship dynamics of two 
identical Mariner hulls. These hulls are superimposed in 
space to allow fcr comparison of the effects contributed to 
rudder modeling differences. In this particular run a step 
and ramp rudder were compared in chapter II. 

Another tenifit of this program is to set up the two 
identical ships reguired for the HAS simulations in chapter 
III. Basically, verification of the models in three degrees 
of freedea is accomplished for the Mariner hull chosen. 

The plots produced in this run are shown in figures 
II-2 and II-2. 
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INSERT TWO /* CARDS HERE 



COMPUTER PROGRAM #2 



This program models a practical rudder response fcr a 
mariner ship type. The rudder limits (stops) are set at ±30 
degrees and the rate of response is limited to ±2 
degrees/seccnd. A scale factor (LUC) is introduced to 
modify the response to match real time of the mariner hull 
chosen. 

Twelve passes thru the program are accomplished to 
conform to different sets of initial conditions and final 
desired rudder conditions. The plots produced in this run 
are shown in figures II-5 and II-6. 
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COMPUTER PROGRAM #3 



This program models a reduced order (first order) gas 
turbine propclsion plant for an input-output relationship. 
The program does not scale the plant to the mariner hull 
used. This was done when introduced into the main 
simulation program first listed as computer program #8. 

The time delay (P) is assisted in initialization by a 
dual feed into the system; one thru the delay itself and one 
directly into SPDIN. The program can be modified to compare 
a family of canes by introducing the following secuence 
into the TERMINAL region: 

INTGER NUMB 

IE (NCUR. EQ. NUMB) CALL ENDEH (NPLOT) 

IF (NCUR. NE. NU-ME) CALL RERUN 
NCUR = NCUR + 1 

where NUMB is the number of curves desired (less than or 
egual to 10) which is set with a PARAM statement. The 
comparison is done on the conditions set in the terminal 
region [i.e. decrement or increment the system gains (eg. G 
= G + 0. C2) ]. 

The plot produced by this run is shewn as part of 
figure 11-10. 
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CCHPOTEE EBOGBAM #4 



This program models a simplified wave simulation 
composed of two superimposed sinusoids (fundamental and 
second harmonic) and a small random wave. The model is 
inherently scaled - to the mariner nondimensional 
characterist ics. Introduction of these waves is accomplished 
in computer program #7. Subroutine DEGEAD is shewn in 
appendix A. 

The sea state force plots in the dimensions of the 
three degrees of freedom produced by this run is shown in 
figures 11-17 thru 11-22. 
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COMPUTES PROGRAN #5 



This picgram uses the mariner hull model first 
introduced in computer program #1 and the control system 
designed in chapter III to simulate the approach phase of 
RAS. The sutroutines and functions that are to be inserted 
from appendix A can also be done in object cede by changing 
the word EORTRAN to OEJECT and placing pre-compiled decks in 
the same locations. In fact, due to the long length of 
subroutine SLOPES , this must be done to be able to run the 
simulation with the DSL default job control language (JCL) 
presently installed at the Naval Postgraduate School IBM 
36C/67 . 

The plots produced by this run are shown in figures 
III-7 thru III-13. By changing the gains and introducing 
the following code, the plots of figures III-14 thru III-19 
are produced: 

CUES = 5 . 0*STEP (8.0)-5.0*STEP (9.0) 
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COMPUTER PROGRAM #6 



This program combines the approach and turn phases of 
computer program #5. The added subroutine is a result of 
simulation requirements to switch between adaptive gains. 

This rcr produced the plots of figures III-22 thru 
III-34. By substituting the initial conditions of table 
III-3, this program produced the plots of figures III-35 
thru IIX-64. 
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COMPUTES EHOGRAM #7 



This program combines the calm sea simulation of 
computer program #6 with the wave simulation of computer 
program 4 4 to simulate the model and control system in a sea 
state. The waves are introduced thru the rudder 
ncndimensicnalized coefficients as shown in chapter II. 

The plots produced are shown in figures III-66 thru 
III-73. figure III-65 was produced with the same program by 
setting fil=1.5. 
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101 FORMAT! 1 THIS RUN IS FOR A PORT SIDE TO APPROACH') 
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102 FORMAT! • THIS RUN IS FOR A STBD SIDE TO APPROACH') 
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INSERT TWO /* CARDS HERE 



COMPOTES PROGRAM #8 



This program incorporates a fifth order polynomial 
curve fit speed control switching function to give optimal 
longitudinal positioning. The scenario is the same that was 
used in the design of the heading control development. The 
low order model of the gas turbine propulsion plant was 
used . 

The plots - produced are shown in figures III-8C thru 
III-83. 
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COMPOTES PHOGBAM #9 



This program introduces a longitudinal position offset 
capabilitj. The method takes the control ship to the 
alongside position until 450 seconds into the run. After 
that time, with the ship steadied, the offset position 
desired (XCFSD) is switched to the desired offset. This 
method negates seme of the transient oscillations which 
cause unstable conditions in the approach phase. A 
secondary change is the use of subroutine SWTCHF instead of 
subroutine 5KTCH developed in the heading control section. 
This new subroutine relaxes the heading velocity feedback 
gain (7F EG) tc allow turn stability in the turn phase. 

The plots produced by this program are shown in figures 
III-84 thru III- 101. 
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COMPUTES PROGRAM #10 



This program incorporates the sea state first 

programmed in computer programs #4 and #7. The WX wave 
force, however, is introduced at the end of the propulsion 
loop- to allow more realistic perturbations. This is the 
final form of the complete heading and speed control 
systems. To run this without a sea state, set (JFMA to 0.0. 

This program produced the plots shewn in figures 
III-103 thru III-105. 
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